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1 .  INTRODUCTION 


This  report  presents  our  experimental  and  theoretical  studies  of 
compression  and  shear  wave  propagation  in  explosives.  The  objective  of 
the  experimental  study  was  to  establish  and  investigate  the  occurrence  of 
shear- induced  reaction  (SIR)  in  an  explosive  subjected  to  one-dimensional 
compression  and  shear  (1DCS)  loading.  The  objective  of  the  theoretical 
study  was  to  formulate  a  theoretical  treatment  for  the  propagation  of 
compression  and  shear  waves  in  an  explosive. 

SIR  in  explosives  was  postulated  in  the  1940s  by  Bridgman, ^  and  the 
results  of  explosive  sensitivity  tests,  such  as  the  drop-weight  and 
fragment  impact  tests,  provide  some  evidence  for  its  occurrence.  Such 
tests  cannot,  however,  be  used  to  establish  the  occurrence  of  SIR, 
because  they  are  not  well  instrumented  and  the  explosives  are  subjected 
to  different  combinations  of  deformations  that  are  not  amenable  to 

O 

quantitative  analysis.  These  limitations  prompted  Gupta  to  propose  the 
use  of  his  1DCS  wave  technique^ ’ ^  to  study  SIR  in  explosives.  In  Gupta's 
1DCS  wave  technique,  the  gas  gun  is  used  to  perform  parallel- inclined 
impact  (PII)  experiments  with  multiple  Lagrange  particle  velocity  gages 
embedded  in  the  target.  In  a  particular  experiment,  the  magnetic  field 
is  oriented  to  allow  the  gages  to  measure  either  longitudinal  particle 
velocity  (u-^)  histories  or  shear  particle  velocity  (U2)  histories.  The 
shear  loading  can  be  varied  by  changing  the  impact  angle.  These  PII 
experiments  remove  the  limitations  inherent  in  the  sensitivity  tests 
because  they  provide  us  with  the  capability  of  (1)  determining  the 
dynamic  shear  properties  of  shock- loaded  material,®  (2)  establishing  the 
occurrence  of  SIR  and  its  role  in  the  initiation  of  an  explosive,  and 
(3)  understanding  the  influence  of  SIR  on  the  flow  produced  in  an 
explosive  by  oblique  impact. 

Our  experimental  approach  to  SIR  can  now  be  introduced  as  an 
application  of  Gupta's  1DCS  wave  technique  to  explosives .  In  this 
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approach,  longitudinal  particle  velocity  histories  are  measured  in  normal 
and  PII  experiments  designed  to  produce  compression  waves  of  the  same 
strength  that  are  nonreactive  in  the  time  frames  of  the  experiments. 

In  this  case,  because  the  loading  conditions  differ  only  by  a  shear 
component,  significant  differences  between  the  sets  of  particle  velocity 
histories  recorded  in  the  normal  and  PII  experiments  can  be  taken  as 
evidence  for  the  onset  and  occurrence  of  SIR  in  the  explosive  subjected 
to  PII. 

The  explosive  Composition  B3  was  chosen  as  the  test  explosive  for 
our  experimental  investigation  of  SIR  because  it  is  readily  available  and 
because  pucks  for  the  gas  gun  experiments  can  readily  be  fabricated  by 
casting  it  around  a  gage  assembly. 

Our  theoretical  approach  to  the  propagation  of  compression  and  shear 
waves  in  an  explosive  is  based  on  an  analysis  of  the  PII  experiment  with 
the  shock- loaded  target  treated  as  an  elastic  solid.  Because  Lagrange 
gages  are  used  in  our  PII  experiment,  we  first  introduced  the  Lagrange 
strain  and  the  first  and  second  Piola-Kirchhof f  (PK)  stresses  as  the 
natural  variables  for  describing  states  attained  in  the  reference 
configuration.  To  provide  a  basis  for  our  investigation,  we  then 
introduced  the  ideal  flow  produced  in  an  elastic  solid  by  PII,  namely,  a 
pure,  inert  compression  wave  (CW)  followed  by  a  pure,  inert  shear  wave 
(SW) .  A  complete  equation  of  state  (EOS)  for  a  reactive,  elastic  solid 
was  then  formulated  and  used  to  derive  the  thermodynamic  equations 
governing  the  adiabatic  flow  produced  in  such  a  solid  by  PII.  Conditions 
for  an  elastic  solid  to  support  the  ideal  flow  considered  earlier  were 
then  derived  and  incorporated  into  a  particular  EOS  for  a  reactive  solid 
so  that  we  could  develop  the  most  fundamental  treatment  of  reactive  CWs 
and  SWs  produced  by  PII.  Flows  with  compression- induced  reaction  (CIR) 
initiated  by  the  CW,  as  well  as  flows  with  SIR  initiated  by  the  SW,  were 
then  considered  to  show  how  the  chemical  energy  influences  the  longitu¬ 
dinal  particle  velocity  and  thereby  provide  a  means  of  interpreting  the 
results  of  the  gas  gun  experiments. 

Section  2  of  this  report  describes  our  experimental  techniques  used 
to  measure  compression  and  shear  response  in  gas  gun  experiments  and  also 
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presents  the  theoretical  background  for  the  PII  experiments.  Section  3 
presents  our  theoretical  study  of  the  reactive  flow  produced  in  a  model 
elastic  solid  subjected  to  PII.  Our  theoretical  study  is  presented 
before  our  experimental  study  because  some  of  its  results  are  used  to 
interpret  the  longitudinal  particle  velocity  histories  recorded  in  the 
gas  gun  experiments.  Section  h  describes  the  experimental  methods  used 
to  fabricate  the  Composition  B3  targets  and  also  presents  the  results 
of  the  gas  gun  experiments.  Our  conclusions  and  recommendations  are 
presented  in  Section  5. 
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2.  THE  PARALLEL -INCLINED  IMPACT  EXPERIMENT 


2.1  EXPERIMENTAL  TECHNIQUES 

This  section  describes  the  techniques  used  to  perform  our 
compression- shear  (CS)  experiments  in  the  gas  gun.  In  these  experiments, 
multiple  Lagrange  gages  are  used  to  record  particle  velocity  histories  in 
a  target  that  is  subjected  to  a  parallel- inclined  impact.  Figure  2-1 
shows  a  schematic  view  of  a  flyer  about  to  strike  a  target  in  a  parallel- 
inclined  impact  experiment.  Upon  impact,  compressive  and  shear  distur¬ 
bances,  parallel  respectively  to  the  directions  X-^  and  X2 ,  are  produced 
at  the  impact  surface  because  the  normal  to  the  flyer  plate  X-^  is  not 
parallel  to  its  direction  of  motion  X^ .  These  compressive  and  shear 
disturbances  propagate  into  the  target  and  back  into  the  flyer  and  will 
be  one -dimensional  with  respect  to  X-^  until  they  are  influenced  by 
lateral  rarefaction  waves.  In  a  normal -  impact  experiment,  when  9=0  and 
X-^  and  X^  coincide,  only  compressive  disturbances  are  produced  at  the 
impact  surfaces.  These  compressive  disturbances  propagate  into  the 
target  and  back  into  the  flyer  and  will  be  one -dimensional  with  respect 
to  X^  until  they  are  influenced  by  lateral  rarefaction  waves. 

The  Lagrange  electromagnetic  particle  velocity  (EPV)  gages  embedded 

in  the  target  are  based  on  the  principle  that  an  electric  conductor  of 
— ♦  — ► 

length  1  moving  with  a  velocity  u  in  a  fixed  magnetic  field  B  will 

generate  a  motional  electromotive  force  E  given  by  the  equation 
— ►  — ¥  — ¥ 

E  =*  1  •  (u  x  B) .  It  follows  from  this  equation  that  we  can  orient  the 
magnetic  field  to  make  either  longitudinal  particle  velocity  (u-^)  or 
shear  particle  velocity  (U2)  measurements. 

A  detailed  account  of  the  experimental  facility  to  produce  and 
measure  properties  of  one-dimensional  compression  and  shear  waves  is 
given  in  Reference  6.  A  discussion  of  theoretical  aspects  of  compression 
and  shear  wave  propagation  is  presented  later  in  this  section  of  the 
report . 
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Figure  2-1 .  Schematic  view  of  experimental  technique  to  produce 
compression  and  shear  waves. 

The  X't  refers  to  the  laboratory  systenv.X,  refers  to  the  coordinate 
system  fixed  on  impact  plates. 


2.1.1  Compression  Response 

To  measure  compression  response  in  a  target,  the  angle  in  Figure  2-1 
is  set  to  zero  in  the  normal  impact  experiments  and  to  15°  or  20°  in  the 
PII  experiments,  and  the  magnetic  field  is  aligned  to  measure  u^(t).  In 
the  assembly  procedure  before  a  gas  gun  experiment,  each  target  is  potted 
in  a  Lucite  target  ring.  Copper  foils  2  mil  thick  are  bonded  to  the 
impact  surface  across  the  target-Lucite  interface,  and  copper  foils  1  mil 
thick  are  bonded  to  the  Lucite  target  ring.  The  2-mil-thick  copper  foils 
are  used  to  trigger  the  oscilloscopes  in  a  gas  gun  experiment,  and  the 
1-mil-thick  copper  foils  are  used  to  measure  impact  alignment. 

The  voltage  signals  from  the  EPV  gages  are  recorded  using  the  usual 
methods  by  analog  and  digital  oscilloscopes.  To  obtain  good  time 
correlations  between  gages,  the  voltages  from  successive  gages  are 
recorded  on  fast  dual-beam  analog  oscilloscopes  (Tektronix  7844) .  This 
procedure  permits  an  accurate  determination  of  wave  velocities  and  an 
improved  analysis  of  the  experimental  data. 

2.1.2  Shear  Response 

To  measure  shear  response  in  a  PII  experiment,  the  angle  9  in 
Figure  2-1  is  set  to  15°  or  20°  and  the  magnetic  field  is  aligned  so  that 
the  voltages  generated  by  the  EPV  gages  have  only  U2  contributions .  Even 
with  good  field  alignment,  shear  particle  velocity  measurements  are  not 
as  accurate  as  longitudinal  particle  velocity  measurements.  The  large 
relative  magnitude  of  the  faster  traveling  longitudinal  wave  coupled  with 
small  errors  in  magnetic  field  alignment  and  impact  tilt  produce  a  small 
perturbation  that  occurs  before  shear  wave  arrival. ^  The  large  relative 
magnitude  of  the  compression  wave  is  a  consequence  of  material  behavior, 
and  merely  increasing  6  does  not  produce  larger  shear  deformation  in  the 
specimen. 

The  target  assembly  for  the  shear  wave  experiments  was  identical  to 
that  for  the  plane  compression  experiments. 
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2.2  THEORETICAL  BACKGROUND 


Our  treatment  of  the  PII  experiment  designed  by  Gupta^’^  is  based  on 
the  assumption  that  the  flow  produced  by  the  impact  is  one -dimensional 
with  respect  to  the  direction  normal  to  the  plane  of  impact.  It  is 
convenient  first  to  introduce  the  strains  and  stresses  that  are  used  in 
our  analysis  of  the  PII  experiment. 

2.2.1  Specification  of  Strain 

Although  we  have  already  used  and  X!^  to  denote  distances  in 
Figure  2-1,  it  is  convenient  at  this  stage  to  introduce  a  new  notation 
and  let  h-^,  x-^,  and  X2  denote,  respectively,  the  Lagrange  coordinates 

and  the  Eulerian  coordinates  normal  and  parallel  to  the  plane  of 
impact.  In  this  case,  to  account  for  the  fact  that  a  particle  may  start 
to  move  with  a  velocity  u-^  in  the  x^  direction  and  a  velocity  U2  in  the 
X2  direction  at  different  times,  we  write  the  equations  for  the  particle 
paths  as 


t 

*1  -  hi  +  /  Ujdt  (1) 

rl 

t 

X2  "  h2  +  u2dt  (2) 

T  2 

so  that  u-^  =*  0  when  t  <  and  U2  -  0  when  t  <  r^.  Remembering  that  the 
flow  is  one -dimensional  and  that  u^  and  U2  are  functions  of  h-^  and  t,  we 
can  differentiate  Equations  (1)  and  (2)  to  obtain  the  following  equations 
for  the  deformation  gradients: 


3x^ 

din 


S 


f ^ 

ah. 


dt 


ax , 


ah. 


(3) 
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3Xr 

Shi 


; 


t  3u2 
Shi 


dt 


dx„ 

_ i 

ah. 


Following  the  notation  used  by  Thurston,^  we  can  then  write  the  matrix  of 
the  deformation  gradients  F  as 


'(dxj/ah^,  o' 
,<ax2/ah2),  iy 


(4) 


and  construct  a  strain  tensor  V  from  the  equation 


V  -  |  (C  -  I) 


(5) 


where  I  denotes  the  unit  matrix,  and  the  right  Cauchy-Green  tensor  C  = 
F+F,  where  F+  is  the  transpose  of  F.  We  can  then  use  the  equations  in 
(3)  to  obtain  the  elements  of  V  describing  the  strain  produced  in  our  PII 
experiments  as 


V 


11 


-  *((?)’  • 


1  (!0z) 

2  \3h1/ 


-  1 


V  -  V 
12  21 


V22  -  ° 


It  is  convenient  at  this  stage  to  introduce  the  symbol  J  for  the 
determinant  of  F, 


(6) 


J  *  det  F 


(7) 


and  present  the  following  equations  for  the  elements  of  C  ^  that  are  used 
later  in  our  discussion  of  stress  in  the  PII  experiment: 
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5h1 

1 

"N 

X 

ro 

-J 

1 

9x^ 

“  J 

^h^ 

3h1 

i  3xi 

I  f _ 

3x2 

J  Uty 

3h2 

1 

3x^ 

j  l-ah^ 

dh2 

1 

dx.. 

r _ = 

3x2 

”  J 

UhxJ 

(8) 


2.2.2  Mass  Balance  Equations 

We  recall  that  the  mass  balance  equation'7  can  be  written  in  terms  of 
the  density  p  or  the  specific  volume  v  =  1/p  as 

pj  -  pqI  v  -  vqJ  (9) 

where  the  subscript  o  denotes  the  initial  state  of  the  material.  It  then 
follows  from  Equations  (4),  (7),  and  (9)  that  the  mass  balance  equations 
take  the  following  forms  for  our  PII  experiment: 

,3x,  dX-  a  5ui 

'feq)  -  V  v  -  vo  feqf  Ji-  v„(5hJ)  (10) 

2.2.3  Specification  of  Stress 

The  first  and  second  Piola-Kirchhof f  (PK)  stress  tensors  are  used  in 
our  treatment  of  material  subjected  to  PII  because  Lagrange  gages  are 
used  in  the  PII  experiments  to  record  particle  velocity  histories  in  the 
reference  configuration. 

The  first  PK  stress  tensor  is  conveniently  introduced  by  trans¬ 
forming  the  equation  for  the  infinitesimal  force  transmitted  in  a  surface 
element  in  the  present  configuration  with  coordinates  (x-^,  X2 ,  x^)  into 
an  equation  in  the  reference  configuration  with  coordinates  (h^,  hj). 
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Q 

To  perform  this  transformation,  we  use  the  exterior  calculus  of  Cartan 
to  write  the  equation  for  a  surface  element  da  with  unit  normal  n  in  the 
current  configuration  as 

da  -  dx^Adx2  +  dx£Adx^  +  dx^dx^  -  nda  (11) 

and  the  equation  for  a  surface  element  dA  with  unit  normal  N  in  the 
initial  configuration  as 

dA  -  dh^Adh^  +  dh^dhj  +  dh^Adh^  e^  -  NdA  (12) 

where  e^,  <  and  e3  are  the  unit  vectors  in  the  current  and  reference 

configurations . 

It  follows  from  Equations  (11)  and  (12)  that  we  can  write 
n^da  -  dx^Adx^ 

n^da  -  dx^Adx^  (13) 

n^da  -  dx^Adx£ 

NxdA  -  dh2Adh3 

N2dA  =  dh^dh  (14) 

N3dA  -  dhXAdh2 

Transfonning  Equation  (13)  into  the  initial  configuration  then  leads  to 
the  equations 

3h-  3h2  3h3 

dx, Adx.  =  J  f- —  dh„Adh-  +  - —  dh.Adh,  +  - —  dh  Adh0] 

1  2  2  3  3x3  3  1  dx^  1  2-> 

3h^  3h2  3h3 

dx2Adx3  -  J  (—  dh2Adh3  +  —  dh^dt^  +  —  dh^dhj  (15) 

3hx 

dx.Adxn  -  J  [- — 

3  1  ^3x2 

and  it  follows  from  Equations  (13)  and  (14)  that  we  can  express  the 
equations  in  (15)  by  the  equation 


3h2  3h3 

dh„Adh,  +  r —  dh-Adh,  +  - — 

2  3  3x2  3  1  3x2 


dhXAdh2] 
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(16) 


n.da  -  J 
J 


ah. 

_ k 

ax. 

j 


N.  dA 
k 


Taking  the  exterior  derivative  of  the  equations  in  (16)  also  leads 
directly  to  the  expression 


Q 

which  is  equivalent  to  the  identity  of  Euler,  Piola,  and  Jacobi. 


(17) 


We  are  now  in  a  position  to  introduce  the  first  PK  tensor  with 
Cartesian  coordinates  Pj^.  We  first  use  the  Cauchy  stress  tensor  with 
Cartesian  components  to  write  the  equation  for  the  infinitesimal 

force  dF  transmitted  in  a  surface  element  da  in  the  current  configuration 
as 


dF  -  a .  .  n .  dae  . 
—  ij  J-  J 


(18) 


We  then  combine  Equations  (15)  and  (16)  to  obtain  the  expression  for  dF 
in  the  reference  configuration  as 


ah. 

dF  -  J  a .  .  ~  ~  N.  dAe . 
—  ij  aXi  k  j 


(19) 


and  set 


ah. 

J  a .  .  ^ 

ij  ax. 


(20) 


It  follows  from  Equations  (18)  and  (19)  that  we  can  compute  the  net 
j -component  of  force  resulting  from  stress  on  a  body  by  integrating 
(7^jn^  over  its  bounding  surface  in  the  present  configuration  or  by 
integrating  over  its  bounding  surface  in  the  reference 

configuration.  Consequently,  in  the  absence  of  body  forces,  the 
equations  of  motion  expressing  the  balance  of  linear  momentum  in  the 
current  configuration, 
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du. 

1 


da .  . 

_ li 

3x . 

J 


(21) 


can  be  transformed  into  the  equations  of  motion  expressing  the  balance  of 
momentum  in  the  reference  configuration. 


3u. 
po  3t 


3P.  . 

ah. 

j 


(22) 


A  consideration  of  the  rate  of  work  W  performed  at  the  bounding 
surface  of  a  body  in  the  current  and  reference  configurations  is  used  to 
define  the  second  PK  stress  tensor.  It  follows  from  our  previous 
discussion  that  W  can  be  expressed  in  the  current  and  reference 
configurations  by  the  equations 


W  -  f  u .  a . . n . da  -  f  U.P..N.  dA 
J  j  ij  l  J  j  jk  k 


(23) 


which  can  be  transformed  by  the  divergence  theorem  into  the  equations 

da . .  au.  ap.,  au. 

_L 

>'i  -*-j  *  *•  •*  ~  j  j** 


w  ,  ,  ui.ii  uv*  . 

"  ‘  f  (uj  ST1  +  "ij  Hpd!VdVd’‘3  -  f  ("j  ShJ-  +  Pjk  ^)<*V«Vdh3 

(24) 


While  Equations  (21)  and  (22)  show  that  the  first  terms  in  these 
integrals  represent  power  of  translation,  the  second  terms  represent 
power  extended  in  stretching  the  material.  The  second  PK  tensor  is 
defined,  to  make  the  equations  for  the  stretching  power  per  unit  volume 
in  the  current  and  reference  configurations  have  the  same  form,  as 
follows.  Because  er—  -  a  j  ^ ,  we  first  write  the  equation  for  the 
stretching  power  per  unit  volume  in  the  current  frame  as 


a .  . 


a .  . d.  . 

1J 


(25) 


where  djj  is  the  stretching  tensor  defined  by  the  equation 


d.  . 
ij 


^  3u^  3u. 

2  tx.  +  3x. 
J  i 


(26) 
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We  next  use  the  equation  dx^Adx2Adx^  =  J  dh^Adt^Adh^  and  the  identity 

•  *7 

relating  d,-  .  and  the  rate  of  strain  tensor  C..  =  2V.  .  , 

&  3-J  ij  ij 


d.  . 
ij 


ah  ah 

— £  —5.  V 
ax.  ax.  rs 
3-  J 


(27) 


to  obtain  the  equation  for  the  stretching  power  per  unit  volume  of  the 
reference  configuration  as 


J a. .d. . 
ij  ij 


ah  ah 

r  s  • 
— -  — -  y 
ax.  ax.  rs 
J 


(28) 


The  second  PK  tensor  is  then  defined  as 


ah  ah 

m  T  r  s 

T  —  J  a.  .  — —  — — 
rs  3x.  ax. 

1  J 


(29) 


to  make  TrsVrs  the  expression  for  the  stretching  power  per  unit  volume  in 
the  reference  configuration. 


It  follows  from  Equations  (20)  and  (29)  that  the  first  and  second  PK 
tensors  are  related  by  the  equations 


3h 

T  -  P.  j-1 
rs  lr  3x. 

J 

and 


(30) 


dx. 

P  =  — -  T 
ir  3h  rs 
s 


(31) 


We  are  now  in  a  position  to  derive  the  relationships  among  ,  P^,  and 
in  our  PII  experiment.  The  equations  relating  the  components  of  the 
Cauchy  and  first  PK  stress  tensors  are  obtained  from  Equations  (20)  and 
(13)  as 


P11  "  all 


dXj 

P12  =  '  an  ah^  +  J<72i 


(32) 
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P21  ’  a12  ~  a21 
d*2 

P22  a12  3h^  +  a22 


The  equations  relating  the  components  of  the  first  and  second  PK  tensors 
are  obtained  from  Equations  (31)  and  (3)  as 


P 


P 


P 


P 


11 


1? 


21 


22 


3x^ 

3h^  T!1 


3x. 

— L  -p 
3h1  21 


3x? 

— -  T  +  T 

ah.  11  12 


3x„ 

-  T  +  T 

3h1  21  22 


(33) 


2.2.4  Equations  of  Motion 

It  follows  from  Equation  (21)  that  we  can  write  the  equations  of 
motion  in  the  current  configuration  for  our  one -dimensional  flow  produced 
by  PII  as 


P 


P 


da 


11 


3x, 


da 


12 


3x, 


(34) 


The  combination  of  Equation  (34)  with  Equations  C9)  and  (33)  then  gives 
the  corresponding  flow  equations  in  the  reference  configuration  [Equation 
(22)]  as 


3o  at 

<3u2 
3o  aF 


a  p 


li 


ah. 


3P 


21 


ah. 


(35) 
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2.2.5  The  Energy  Equation  and  Equations  of  State 


To  complete  the  balance  equations  for  our  PI1  experiment,  we  assume 
that  the  flow  in  material  subjected  to  PII  is  adiabatic.  In  this  case, 
it  follows  from  the  first  law  of  thermodynamics,  with  W  specified  by 
equation  (23),  that  we  can  write  the  equations  for  the  time  rate  of 
change  of  the  specific  internal  energy  e  in  the  current  and  reference 
configurations  as 


P 


de 

dt 


a .  .d. . 


and 


de 
p o  dt 


3u. 
P  — - 

ij 

J 


de 
?o  3 1 


-  T. .V. . 
ij 


(36) 


(37) 

(38) 


It  then  follows  from  Equations  v37)  and  (38)  that  we  can  write  the  energy 
equation  for  our  PII  experiment  as 


de 
Po  dt 


du^  3u2 

P11  ih^  '  p21  3h^ 


(39) 


or 


de 
1  o  dt 


-TV  -TV  -TV 
11  11  12  12  21  21 


(40) 


To  provide  a  more  complete  description  of  the  flow,  we  assume  that 
the  material  subjected  to  PII  is  an  inert  elastic  solid  and  denote  its 
specific  entropy  and  temperature  by  s  and  T.  The  differential  forms  of 
the  e  -  e(s ,  vQ3x^/3hj )  and  e  =  e(s,vQV^j)  equations  of  state  (EOS)'*'®’'7 
can  then  be  written  as 


de  =  Tds 

-  v  p. ,d(3x./3h. ) 
oKij  x'  j 

(41) 

de  =  Tds 

-  v  t .  .  dV .  . 

(42) 

o  1J 
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where  P^j  =  *  de/d (vQ3x^/3hj )  and  t  —  =  -  3e/3(v0V^j)  are  the  thermo¬ 
dynamic  stresses  corresponding  to  the  first  and  second  PK  stresses. 
Equations  (37)  and  (41)  and  Equations  (38)  and  (42)  can  then  be  combined 
to  give  the  following  equations  for  the  entropy  production  in  adiabatic, 
nonreactive  flow. 

3  3u. 

T3l  '  *  VPij  -  Pij>  3b1 

J 

T  fz  -  -  v  (T. .  -  t. .)  V. .  (44) 

3t  ov  ij  ij'  ij 

Equations  (43)  and  (44)  are  expressions  of  the  second  law  of  thermo¬ 
dynamics  for  the  adiabatic  flow  of  an  inert  elastic  solid.  To  simplify 
our  treatment  of  the  adiabatic  flow  in  a  reactive  elastic  solid,  we  will 
assume  that  the  entropy  production  is  dominated  by  the  reaction  and 


accordingly  set  Pi 


p^j  and  -  t^j  in  regions  free  from  shocks. 


2.3  IDEAL  INERT  FLOW  PRODUCED  BY  PARALLEL- INCLINED  IMPACT 

Our  analysis  of  the  ideal  flow  produced  in  a  PII  experiment  is  based 
on  the  assumption  that  inert  material  subjected  to  PII  can  support  two 
waves,  a  compression  wave  (CW)  in  which  d^/dt  -  0  and  a  subsequent  shear 
wave  (SW)  in  which  3u^/3t  -  0.  The  conditions  required  to  obtain  this 
type  of  flow  are  formulated  later  in  this  report. 

2.3.1  Flow  in  the  Pure  Compression  Wave 

When  3u2/3t  -  0  in  the  CW,  we  use  the  identities 


t  3ur 


U2  “  ^  3t 


3u„  t  3^u„ 

rr r  ~  f  dt 

3hn  J  dhndt 

1  'i  1 
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to  show  that  U2  -  0  and  di^/dh^  =0.  It  then  follows  from  Equations  (3) 
(6),  (33),  (35),  and  the  equation 


h  dP  .  h1  3u. 

P2i  ~J  sit-  dhi  '  '  J  a -  dhi 

h. ..  1  h. 

lf  if 


that  the  flow  in  the  CW  also  satisfies  the  equations 
3x„ 


11 

2 

3h 

=  V 

=  0 

12 

21 

21 

=  T 

12 

=  0 

We  will  now  perform  a  Lagrange  analysis  with  the  mass  balance  and 
momentum  equations , 

a  <9u, 

3v  _  _ 1 

dt  Vo  dh^ 

3U1  3P11 


to  determine  stress  -  strain  conditions  for  a  CW  to  form  a  shock  and 

propagate  at  a  constant  velocity.  Extending  the  notation  introduced  by 
11  12 

Fowles ,  ’  we  set  (3h^/3t)u^  -  and  rewrite  Equations  (52)  and  (53) 

in  integral  form  as 


v  -  v  -  v  I  — 
o  o  J  3h 


t  3u. 

f 


-  v  f  1  C_1du 
o  11 


h,  3u 


r  1  1  ,,  r 1  r  a 

“o  .  J  ar  dhi  -  "o  J  cidui 
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Differentiating  Equation  (54)  with  respect  to  hq  at  constant  Uq  and 
Equation  (55)  with  respect  to  t  at  constant  uq  then  leads  to  the 
equations 


(56) 


(57) 


which  show  that  v  -  v(uq)  and  Pq^  -  ^ll^ul^  when  Cq  -  Cq(uq)  . 

We  now  consider  such  a  simple  wave  with  Cq  =■  Cq(uq)  .  We  rewrite 
Equations  (52)  and  (53)  as 


du. 


d(3x^/3h^) 
dP, 


-  C, 


11 


du. 


PoCl 


(58) 

(59) 


and  then  combine  Eolations  (58)  and  (59)  to  obtain  the  equation 
dP, 


'll 


d(3x1/3h1) 


“A 


(60) 


which  relates  the  slope  of  the  first  PK  stress  -  strain  curve  in  the 
reference  configuration  to  the  square  of  the  propagation  velocity  Cq.  It 
follows  from  Equation  (60)  that  the  slope  of  the  stress  (Pqq)  -  strain 
(<3xq/dhq)  curve  is  negative.  Differentiation  of  equation  (60)  with 
respect  to  Uq  then  gives  the  equation 


dC] 

du. 


d2P 


11 


2cJ  d(3x1/8h1)2 


(61) 


which  can  be  used  to  derive  conditions  for  a  simple  CW  to  steepen, 
flatten,  or  remain  unchanged  as  it  propagates.  We  now  note  that  dCq/duq 
>  0  is  the  condition  for  a  simple  CW  to  steepen,  dCq/duq  <  0  the 
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condition  for  a  simple  CW  to  flatten,  and  dC^/du-^  =  0  is  the  condition 

for  a  simple  CW  to  propagate  unchanged.  It  then  follows  from  Equation 

2  2 

(61)  that  a  simple  CW  will  steepen  to  form  a  shock  if  d  PixA^Sx^/Sh^) 

>  0,  will  flatten  and  spread  out  if  d^P^/dOx^/Sh-^) 2  <  0,  and  will 

2  2 

propagate  unchanged  if  d  P^/d^x^/Sh-^)  =  0. 

We  now  use  Equations  (33)  and  (49)  to  formulate  these  conditions  in 
terms  of  the  relationship  between  the  second  PK  stress  and  the  V-q  strain 
as  follows.  Differentiating  Equations  (33)  and  (49)  leads  to  the 
equations 


dP 


11 


3x,  2  dT 


d(3x1/5h  ) 


T  + 
11 


(ar) 


dV 


11  -  „  r2 

P  '-'-i 

n  0 1 


and 


d2P 


11 


3x, 


d(3x1/3h  )‘ 


fer) 


dT 


11 


3h J  dV 


11 


3x1 

fe1) 


2 

3  dT 


11 


aV  W2 


(62) 


(63) 


It  follows  from  Equation  (62),  because  T^  >  0  and  dT^/dV^  <  0, 
that  the  (T-q  -  V-^)  stress-strain  curve  must  satisfy  the  condition 
IdTu/dVnl  >  Tn(3hi/3xi)^ .  It  also  follows  from  Equation  (63)  that 
conditions  for  shock  formation  are  more  complicated  when  they  are 

expressed  in  terms  of  the  (T-^  -  V^)  curve  rather  than  the  (P^  - 

2  2 

3x-^/3h^)  curve  because  the  condition  d  P^/d(3x^/3h^)  >  0  ensures  shock 

2  2 

formation  in  a  simple  CW  but  the  condition  d  Ti]_/dV-^  >  0  does  not. 

Steady-State  Flow.  We  now  consider  the  situation  when  the  CW 
attains  a  steady- state  and  propagates  unchanged  with  a  velocity  D.  In 
this  case,  because  dC-^/du^  =  0 ,  we  can  set  C^  =  D  and  integrate  Equations 
(58),  (59),  and  (60)  to  obtain  the  Rankine -Hugoniot  jump  conditions 


(64) 

(65) 

(66) 
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We  can  then  use  the  following  equation  obtained  from  Equations  (49)  and 
(33) 


pn  -  Ki + 1)1/2  Tn  (67) 

to  express  Equations  (65)  and  (66)  in  terms  of  the  second  PK  stress 
tensor  as 


11  -  poDul 

(2Vu  + 1) 

-1/2 

(68) 

11  - 

(2vn  +  0 

-  l) 

(69) 

We  can  then  use  Equation  (66)  to  integrate  the  energy  equation 

3u, 

de  _ 1 

Po  3t  “  *11  6h, 


(70) 


to  obtain  the  Hugoniot  relation 


3x, 


6  -  eo  -  2  VllO  -  3h^) 

which  is  readily  transformed  into  the  equation 

!C  -  T  Tn((2vu  *  x)1/2  •  On  +  l)) 


(71) 


e  -  e 


(72) 


Thermodynamic  Considerations.  We  are  now  in  a  position  to  use  the 
e  =  e(s , v03x^/3hj )  and  e  =  e(s,v0V^j)  EOS  to  consider  thermodynamics 
aspects  of  the  flow.  At  this  stage,  for  convenience  in  writing  partial 
derivatives,  we  set  3x^/3h^  =  x^  and  dy.^/dh^  =  x£ .  For  this  purpose,  it 
is  convenient  to  discuss  the  isentropic,  simple  CW.  In  this  case,  we  can 
set  -  p^j  and  T^j  -  t—  in  Equations  (43)  and  (44) ,  rewrite  Equations 
(60)  and  (61)  as 


3p 


11 


^3(3x1/ah 


-) 

l  s,x' 


=  *  %C1 


(73) 


21 


dC1  v 
1  o 


I  "  pu  \ 

i  2c2  VafavahJ2/ 


and  Equations  (62)  and  (63)  as 


t, ,  + 


(2vii + »  fe»r) 


11  s  V  V 
’12*  21 


3(5x1/ah1)‘ 


wo  3t- 1 

3(2V  ♦  1)1/2  [-H) 


11  S'V12'V21 


+  (2V1L  +  1) 


11  s  V  V 

’12’  2! 
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where  (3t^/3V^) s  ^  and  (3  t^/3V^) g  ^  are  respectively  called  the 
second-order  and  thircl-order  isentropic  elastic  stiffness  coefficients. 
Equations  (10)  and  (35),  written  for  isentropic  flow,  can  then  be 
combined  to  show  that  the  slopes  of  the  characteristics  in  the 
(t  -  h-^)  plane  are  given  by  the  equation 


-  +  -v 


o  3(3x 


pn  7 

1/3V  s.J 


and  it  then  follows  from  Equation  (73)  that,  in  a  simple  isentropic  CW, 
the  lines  of  constant  particle  velocity  in  the  (t  -  h^)  plane  are  also 
the  forward  facing  characteristics.  Equation  (74)  for  the  change  in  the 
slopes  of  the  lines  of  constant  particle  velocity  with  particle  velocity 
is  thus  equivalent  to  the  equation  given  by  Courant  and  Friedrichs  for 
the  change  in  the  slopes  of  the  characteristics  with  particle  velocity 
for  a  fluid.  In  this  light,  it  is  clear  that  can  be  calculated  in  our 
simple  isentropic  CW  with  3x2/3h^  -  0  from  the  e  =  e(p-j_^,  v^.v^xj)  EOS 
as  follows.  Equating  the  expressions  for  de  obtained  from  the  e  = 
e(s  ,vQx^  ,v  x£)  and  e  -  e(Pll  -  vQxi  ’v0xp  E0S  leads  to  the  equation 
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Tds  -  pud(vo(3V3h1)  -  [£-) 


dp 


11  x^,x£ 
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3e 


l3(v  3x1/3h1  , 
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y)  d(v0(axi/3M 


(78) 


which  reduces  to  the  following  equation  when  ds  -  0: 


r — flii — ) 

la<3xi/ahi)Js> 
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pll’x2 
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(79) 


When  the  e  -  e(p-^^ ,  vQx^ ,  v^x^)  EOS  is  known,  Equation  (79)  can  be 
integrated  to  obtain  the  (p^  -  3x^/3h^)  expression  for  an  isentrope ,  and 
then  Equations  (73)  and  (74)  can  be  used  to  calculate  values  of  C-^  and 
the  corresponding  values  of  dC^/du^. 

We  finally  consider  equations  (73,  (75),  (76),  and  (79).  First,  we 
note  that  Equation  (79)  can  be  transformed  into  Equation  (75)  by  making 
use  of  the  equations 
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(81) 


(82) 


We  next  note,  from  Equations  (73)  and  (75),  that  the  initial  slopes  of 
the  (P]_i  -  3x^/3h^)  and  (t-^  -  V^)  isentropes  are  the  same  and  define 
the  propagation  velocity  of  the  foot  of  a  simple  isentropic  CW.  Finally, 
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we  note  from  Equation  (76)  that  model  elastic  solids  that  do  not  include 
third-order  isentropic  elastic  stiffness  coefficients  in  the  formulation 
of  their  e  -  e(s,v0V^j)  EOS  will  not  support  the  propagation  of  a 
compressive  shock  because  (3  t^/dV^)  g  ^21  < 

2.3.2  Flow  in  the  Pure  Shear  Wave 

Our  treatment  of  the  inert  pure  SW  parallels  our  treatment  of  the  CW 
and  is  based  on  the  assumption  that  the  CW  has  attained  a  steady  state 

A  A  A  A  A  A  A  A 

with  the  peak  values  v  =  VQ(dx^/dh^) ,  u^,  =  p^,  =  t^,  and  e. 

When  dU]_/3t  «  0  in  the  SW,  we  use  the  identities 


a  t  du.. 

U1  -  U1  “  / 


2 

du..  t  d  u. 

— —  -  f  - ±- 

dh..  J  dtudt 

1  r2  1 


to  show  that  u-^  =■  u^  and  du^/dh^  =  0.  It  then  follows  from  Equations 
(3),  (6),  (33),  and  (35)  and  the  equation 


P  -  P 
11  11 


h1  dP..  ..  h..  du 

J  5hT  dhi  '  '  %  J  aT  dhi 
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24 


T11  -  T11  ■  hi  -  'll 


(^r 

\3h  / 


We  now  assume  that  the  SW  is  a  simple  wave  and  write  the  equations 


a  3x2 


3t  3h^  3h^ 


in  terms  of  (Sh^/St)^  -  -  C2(u2>  as 


d(3x2/ah1) 


Following  our  previous  discussion  of  the  CW,  we  can  then  write 
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.  _  .  r> 

d(3x2/3h1) 

-  P0C 

^2m\_ 

d2p21 

dU2  2C2 

d(3x2/5h 

It  then  follows  from  Equations  (95)  and  (96)  that  the  slope  of  the  (P21 

dx2/3h-L)  curve  must  be  negative  and  that  the  SW  will  form  a  shock  if 

d^P2^/d(3x2/3h-^) ^  >  C,  flatten  and  spread  if  d^P22yd(3x2/3h^) ^  <  0,  and 

2  2 

propagate  unchanged  if  d  P2^/d(3x2/3h^)  —  0. 

We  now  use  Equation  (88)  and  the  equation  from  (33) 
dx.  * 

P  =  — -  T  +  T  (97) 

21  ahx  11  12  k 
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to  express  Equations  (95)  and  (96)  in  terms  of  the  components  of  the 
second  PK  stress  tensor.  Differentiating  Equation  (97)  twice  with 
respect  to  (3x2/3h^)  gives  the  equations 


d?21  *  l  ^12  _  2 

d(3x2/ah1)  'll  +  2  dV12  '  P0U2 


d2f21  _  1  ^12 
d(ax2/ah1)2  4  dv22 


We  now  consider  Equations  (96),  (97),  (98),  and  (99).  Because 
(3x2/3h-^)  <0,  it  follows  from  Equation  (97)  that  the  condition  T^2  > 

-  T.^  (3x2/3h^)  is  satisfied  in  the  SW.  Because  >0,  it  follows  from 
Equation  (98)  that  dT^2/dV^2  <  0  and  also  that  the  condition  -  dT-^/dV-^ 

>  2T^  is  satisfied  in  the  SW.  It  follows  as  a  consequence  of  the  latter 
condition  that  the  shear  strain  will  be  located  close  to  the  impact 

A 

surface  if  -  dT-^/dV-^  ~  2T^ .  It  follows  from  Equations  (96)  and  (99) 
that  the  tendency  for  a  SW  to  steepen,  flatten,  or  propagate  unchanged  is 
determined  solely  by  the  sign  of  the  second  derivative  of  the  (T-^  -  V^2 ) 
stress -strain  curve.  It  is  clear  that  the  conditions  for  steepening, 
flattening,  and  propagating  unchanged  are,  respectively,  d2Tj_2/dV22  >  0, 
d2T12/dV12  <  °>  and  d2T12/dV^2  -  0. 

Steady  State  Flow.  We  will  now  derive  the  equations  for  a  shear 
shock,  although  a  previous  study**  leads  us  to  expect  that  solids  satisfy 
the  condition  d2T^2/dV2^2  <  0-  When  d^T-^/dV^  -  0 ,  we  can  set  C2  -  Ds 
and  integrate  Equation  (98)  to  obtain  the  equations 


-  "o  (Sj)  “2 

-  vvi2  (d2  +  °2  P  •  0) 


(100) 


(101) 


We  can  then  use  Equation  (100)  to  integrate  the  energy  equation 


a  _ 

oe  _ / 

po  at  “  21  ah. 


(102) 
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to  obtain  the  equation 


A 


e  -  e 


which  can  be  rewritten  as 


(103) 
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12  11 
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Thermodynamic  Considerations.  Following  our  previous  discussion  of 
the  CW,  we  will  consider  an  isentropic,  simple  SW.  Then  we  can  set  P—  = 
Pij ,  T^j  -  t^j  and  rewrite  Equations  (95)  and  (96)  and  (97)  as 


3<3V-,hi)s,xi  ‘ ' 

dC2  Vo  /  3\l 

du2  2C2  \3(3x2/ah1 


P21 


hi  +  C12 


(105) 
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It  then  follows  by  differentiating  Equation  (107)  and  taking  account  of 
Equation  (87)  that  Equations  (98)  and  (99)  can  be  rewritten  as 
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Now  because  V12  <  0,  (at12/3V11)s f y^  # <  0,  and  (3t12/3V12)s ^  f  < 
0,  it  follows  from  Equation  (108)  that  the  condition  -  (St-^/dV-^) 


>  t^  +  2V^2 s  y  y  is  satisfied  in  a  simple  isentropic  SW. 
Comparing  Equations  (109)  and  (77)  shows  that  the  expression  for  the 
second  derivative  of  p2]_  with  respect  to  (3x2/3h^)  is  more  complicated 
than  the  expression  for  the  second  derivative  of  p-^  with  respect  to 
(3x^/3h-^)  .  This  difference  obviously  arises  because  is  only  a 

function  of  (3xj/3h^)  in  the  CW  but  V-q  and  V-^  are  both  functions  of 
(3x2/3h-^)  in  the  SW.  It  also  follows  from  Equation  (109)  that  model 
elastic  solids  that  do  not  include  third-order,  isentropic,  elastic 
stiffness  coefficients  in  the  formulation  of  their  e  =  e(s,vQV^j)  EOS 
will  support  a  spreading  SW  if  (St^/dV^) s  y  y  ?*  0  but  a  steady- 
state  SW  if  (3t^2/dvii)  -  0. 


The  combination  of  Equation  (92)  ,  written  for  isentropic  flow, 
with  Equation  (91)  gives  the  following  equation  for  the  slopes  of  the 
characteristics  in  the  (t  -  h-^)  plane, 


dt 


+ 


3p21 

Vo  3 (3x2/3h1) s 


(HO) 


Comparing  Equations  (105)  and  (110)  then  shows  that,  in  a  simple 
isentropic  SW,  the  lines  of  constant  particle  velocity  in  the  (t  -  h-^) 
plane  are  also  the  forward  facing  characteristics.  From  an  argument 
similar  to  that  presented  earlier,  it  then  follows  that  the  differential 
equation  for  the  (p2^  '  dx2/dh^)  isentropes  can  be  written  in  terms  of 
the  e  =  e(p2]_,v03x-^/3h-]_,v03x2/3h^)  EOS  as 


3p21 

3(3x2/3h1) 


s,x^ 


V 


P21  +  (3e/3vo3x2/3h1) 


P21’X1 


01  <3e/3«,2i)*p*; 
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Equation  (111)  can  then  be  transformed  into  Equation  (108)  by  making  use 
of  Equation  (105)  and  the  thermodynamic  identities 
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(112) 


(113) 


(114) 


(115) 


3.  THEORETICAL  STUDY  OF  REACTIVE  FLOWS  PRODUCED  IN  AN  ELASTIC 
SOLID  BY  NORMAL  AND  PARALLEL- INCLINED  IMPACT 


Up  to  this  point  in  our  theoretical  treatment,  we  have  assumed  that 
the  flow  produced  in  material  subject  to  normal  or  parallel- inclined 
impact  is  nonreactive.  To  remove  this  restriction,  we  need  to  formulate 
an  EOS  for  a  reactive  mixture  and  then  combine  this  EOS  with  the  energy 
equation  to  determine  how  the  reaction  influences  the  flow. 


3.1  EQUATION  OF  STATE  FOR  A  REACTIVE  EXPLOSIVE  MIXTURE 

The  formulation  of  an  EOS  for  a  reactive  mixture  of  an  explosive  and 
its  products  is  a  complicated  and  controversial  problem.  We  here 
consider  the  simplest  case  by  assuming  that  the  explosive  undergoes  a 
global  decomposition  to  produce  solid  reaction  products  and  a  reactive 
mixture  that  is  in  thermal  and  mechanical  equilibrium  and  obeys  the 
following  mixture  rules  for  the  specific  internal  energy,  the  specific 
entropy,  and  the  strain,  respectively: 


e  =*  (1  -  A)e  +  Ae 
x  p 


(116) 


s  -  (1  - 

A)s  +  As 
x  p 

(117) 
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.  .  x  l  ,  ,  p  1 

(1  •  1>v0  ah.  +  Xvo  8h. 

(118) 

In  these  equations,  we  use  the  subscripts  and  superscripts  x  and  p  to 
denote  the  explosive  and  its  products,  and  we  use  A  to  denote  the  mass 
fraction  of  explosive  that  has  decomposed.  We  now  let  Q  denote  the 
standard  heat  of  decomposition  of  the  explosive,  rewrite  the  specific 
energies  of  the  explosive  and  its  products  as 
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•Vo^xj/ahj)31) 
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=  -Q  +  e°  +  e  fs  ,vP(3x./3h.  )P1 
x  p'-p  o  1  y  ) 


(120) 


and  combine  Equations  (116),  (119),  and  (120)  to  obtain  the  equation 


e  =  e  -  AQ  +  (1  -  A)e  +  Ae 
x  x  p 


(121) 


We  next  combine  the  differential  forms  obtained  from  Equations  (121)  , 
(117)  ,  and  (118)  with  Equation  (41)  written  for  the  explosive  and  its 

X  T) 

products  and  the  equilibrium  conditions,  -  T  and  p..  =  p..,  to  obtain 

x  P  !J  ij 

the  following  differential  form  expressing  e  as  a  function  of  s, 

3x./3 h. ,  and  A : 
x'  j 


de  =■  Tds  -  v  p  .  . d(3x. /3h . )  -  AdA 
o  ij  i  j 


(122) 


where 

dx.  x  _  dx.  p 

A  “  AQ  +  [®x  Tsx  +  p  -  p  -  Ts  +  p  ]]  (123) 

J  J  J 

At  this  stage,  we  have  shown  that  a  reactive  explosive  mixture  obeying 
our  assumptions  can  be  described  by  an  e  -  e(s , vQ3x^/3hj , A)  EOS.  We  must 
now  combine  Equation  (122)  with  the  flow  equations  to  obtain  the 
corresponding  thermodynamic  equations  governing  reactive  flow. 

3.2  THERMODYNAMICS  OF  THE  REACTIVE  FLOW  OF  AN  ELASTIC  SOLID 

Our  derivation  of  the  thermodynamic  equations  governing  the 
adiabatic  reactive  flow  of  an  elastic  solid  is  based  on  the  ZND 
assumption^  that  the  entropy  production  in  regions  free  from  shocks  is 
associated  solely  with  the  chemical  reaction.  In  this  case,  we  can  set 
Pij  ”  P i j  >  as  discussed  previously,  and  combine  Equation  (39)  with 
Equation  (122)  written  as 


3e  3s  ^1  3u2  3  A 

3t  “  3t  VoPll  3hx  ‘  Vop21  3h  ~  3t 


(124) 
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to  obtain  the  basic  equation  governing  our  reactive  flow 
3s  A  3A  r3s~>  3  A 

at  ”  t  at  UaJ  ,  ,  at  { 

e,x|,x£ 

The  combination  of  equation  (125)  with  the  s  -  s(Pll  >v0x-[  .v^ ,  A)  and 
s  —  s (P21 > V0X| . V0X2 > A)  EOS  then  leads  to  the  equations 
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Combining  Equation  (126)  with  the  T  =  T(Pn » VQX^ ,  vqx2 ,  A)  EOS  and  Equation 
(127)  with  the  T  =  T(p2^ , v^x^ , v^x^ , A)  EOS  leads  to  the  equations 
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where  the  specific  heats  C  ,  and  C  ,  are  defined  by  the 
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and  satisfy  equations  of  the  form 
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The  comparison  of  Equations  (126)  and  (127)  with  Equations  (130)  and 
(131)  demonstrates  an  important  feature  of  shock- induced  reactive  flow. 
Because  we  expect  the  ratios  of  specific  heats  in  Equations  (130)  and 
(131)  to  be  approximately  one  and  the  derivatives  (3p  /3A)  ,  , 

11  L  t X  ^  , X  ^ 

and  (3p91/3A)  ,  ,  to  be  small  compared  to  the  derivatives 

4  A  —  > x^ 1 ^0 

(3p11/3A)  ,  ,  ana  (3p91/3A)  ,  ,  ,  the  temperature  can  increase  along 

ii  e , x^ , X2  z  i  e ,x^ ,X£ 

a  particle  path  as  the  stresses  decrease.  Consequently,  for  a  predomi¬ 
nantly  temperature  dependent  reaction  rate,  the  rate  of  reaction  can 
increase  along  a  particle  path  as  the  stresses  decrease,  and  the  shock 
initiation  of  detonation  can  occur  in  a  wave  satisfying  the  condition 
3pil/3t  <  0. 
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We  next  consider  the  e  -  e(T,v  x'  v  x'  A)  EOS  and  write  its 

o  1  o  l 

differential  form  as 
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To  interpret  the  derivatives  (3e/3v  x' )_  ,  and  (3e/3v  x')T  we 

Ol  i , X2  f  A  O  2  I , X^ , A 

rewrite  Equation  (136)  for  a  nonreactive  fluid  as 

™  .  f  3e  1  _ 


de  -  C  dT  +  dv  (137) 
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O  i  T 

and  use  the  first  law  of  thermodynamics  for  adiabatic  flow  written  as 

de  -  -  5W  (138) 


to  obtain  the  equation 
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(139) 


We  now  recall  that  (3e/3v  x' )„  is  associated  with  forces  of  interaction 

oil 

among  molecules  and  use  Equation  (139)  to  show  that  the  work  5W  performed 
on  a  fluid  is  partitioned  between  thermal  energy  and  repulsive  energy 
when  (3e/3vQX^)^  <0.  We  also  recall  that  for  a  gas  with  (Se/Sv^x^)^  =  0, 
all  the  work  is  converted  into  thermal  energy,  but  for  a  liquid  with 
(3e/3vQxp^,  <0,  a  large  fraction  of  the  work  can  be  converted  into  the 
repulsive  energy  arising  from  compression.  In  Equation  (136),  we  thus 
associate  (3e/3v  x' )  and  (3e/3v  x’)  ,  .respectively,  with  the 

O  1  I , X2  j  A  O  2  X.  , A 

repulsive  forces  arising  from  compression  ana,  shear  and  on  physical 
grounds,  assume  that  (3e/3v  x' )  ,  /(3e/3v  x')  ,  >  1. 

O  1  X  1  X 2  )  A  O  2  i  ,  X^  ,  A 

We  finally  combine  Equation  (136)  with  Equation  (39)  to  obtain  the 
equation 
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3.3  MODEL  REACTIVE  SOLID  FOR  THE  GAS  GUN  EXPERIMENTS 


To  include  chemical  reaction  in  our  treatment  of  the  PII  experiment, 
we  adopt  the  classical  ZND  model  of  detonation^  and  assume  that  the 
front  of  the  wave  produced  in  such  an  experiment  is  a  nonreactive  shock 
discontinuity.  Because  we  are  interested  primarily  in  the  dependence  of 
the  longitudinal  particle  velocity  on  the  chemical  reaction,  it  is  conve¬ 
nient  to  restrict  our  attention  to  an  equation  of  state  for  the  reactive 
mixture  that  provides  the  most  basic  treatment  of  this  dependence. 
Accordingly,  we  chose  an  equation  of  state  that  is  compatible  with  the 
ideal  inert  flow  discussed  in  Section  2.3  above.  Properties  that  must  be 
exhibited  by  this  equation  of  state  are  derived  as  follows. 

Conditions  for  the  shear  particle  velocity  in  the  reactive  shock 
discontinuity  at  the  wave  front  to  remain  zero  can  be  inferred  from  the 
following  equations,  obtained  respectively  by  setting  dX/dt  =  0  in 
Equation  (127)  and  differentiating  Equation  (47): 


It  follows  from  Equations  (127)'  and  (141)  that  the  assumption  that 
dn^/dt  -  0  in  a  leading  inert  isentropic  CW  is  only  truly  valid  for 
materials  that  satisfy  the  condition  (dp91/dx')  ,  -  0.  Now  because 

Z,  X.  J.  S  y  X  2 

du^/dh^  -  0  when  du^/dt  =0,  it  follows  from  Equation  (97)  written  for 
the  thermodynamic  stresses  that 
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We  thus  conclude  that  material  subjected  to  PII  will  support  a 
nonreactive  CW  with  no  shear  velocity  component  when  at^/dV^  is 
proportional  to  V^- 
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Conditions  for  the  condition  3u^/3t  -  0  to  be  satisfied  in  an 
ensuing,  inert  isentropic  SW  can  be  inferred  from  the  following 
equations,  obtained  respectively  by  setting  3A/3t  -  0  in  Equation  (126) 
and  differentiating  Equation  (85) : 


3pn 

at 


S  , 


apn 

3t 


dt 


(126)' 


(143) 


It  follows  from  Equations  (126)'  and  (143)  that  the  assumption  that 
3u^/3t  =  0  in  an  ensuing  inert  isentropic  SW  is  only  truly  valid  for 
materials  that  satisfy  the  condition  (3p1-/3x‘)  ,  -  0.  The  equation 
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obtained  by  differentiating  the  equation  relating  p-Q  and  t^  [Equation 

(33)],  shows  that  the  condition  (Sp^/Sx^)  s  x<  -  0  is  much  more 

restrictive  than  the  condition  (3p91/3x')  ,  i  0. 

Z  J-  L  S f ^2 

It  follows  from  the  above  discussion  that  the  ideal  inert  flow 
described  in  Section  2.3  can  be  modeled  by  using  an  EOS  satisfying  the 
following  conditions: 
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A  simple  model  EOS  that  has  the  initial  hydrostatic  stress  pQ  -  0  and 
satisfies  these  conditions  is  conveniently  written  in  terms  of  x£  and  the 
displacement  gradient,  u'  -  3u/3h^  -  x^  -  1,  as  follows: 
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It  follows  from  the  definition  of  u'  that  we  can  set  (d/du')  -  (5/<3x^) 
when  we  differentiate  Equation  (146).  Differentiating  Equation  (146) 
partially  with  respect  to  A,  s,  u' ,  and  leads  to  the  following 
equations  for  A,  T,  p^]_,  and  p2]_: 
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It  is  clear  from  these  equations  that  we  are  neglecting  the  difference 

between  A  and  Q,  the  initial  temperature,  and  the  derivatives 

(dp01/3A)  ,  ,  and  (3p91/3A)  ,  Some  Maxwell  relationships  that 

“  -*»  S  j  f  4  1  s  ,  f  Xrt 

we  will  use  later  are  obtained  by  differentiating  Equations  (148) ,  (149) , 
and  (150)  as 
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It  is  also  clear  from  Equations  (149)  and  (150)  that  our  model  EOS 
satisfies  the  last  two  conditions  expressed  in  equation  (145).  For 
consistency  with  the  theory  of  elasticity,  we  set  A-^(s)  =  6  +  Ip,  A2(s)  -  n 
and  note  that  5  is  used  in  the  expression  for  A-^(s)  because  A  has  already 
been  used  for  the  reaction  coordinate.  Equations  (149)  and  (150)  then 
give  the  following  governing  equations  for  our  ideal  inert  flow  in  a  PII 
experiment : 
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After  combining  Equations  (73),  (74),  (155),  and  (157)  and  Equations 
(105)  ,  (106) ,  (156)  ,  and  (157) ,  we  conclude  that  in  the  absence  of  a 
chemical  reaction,  the  isentropic  discontinuity  at  the  front  of  the  pure 
CW  propagates  at  the  velocity  C-^  -  [  vQ(5  +  2/i)]^2,  and  the  isentropic 
discontinuity  at  the  front  of  the  ensuing  pure  SW  propagates  at  the 
velocity  C2  -  . 

We  are  now  in  a  position  to  consider  our  model  for  the  reactive  flow 
produced  in  a  PII  experiment.  The  combination  of  the  time  derivative  of 
equation  (146)  with  the  first  law  of  thermodynamics  expressed  by  Equation 
(39)  leads  to  the  following  equation  for  the  entropy  production 
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which  shows  that  the  only  irreversible  process  being  considered  is  the 
chemical  reaction.  The  other  equations  governing  our  reactive  flow  can 
readily  be  derived  as  follows. 

We  combine  the  equations  obtained  by  differentiating  Equations  (149) 
and  (150)  with  Equations  (148)  and  (158)  to  obtain  the  time  derivatives 
of  and  p2]_  as 
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It  then  follows  from  Equations  (126)  and  (127)  that  for  our  EOS 
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To  obtain  the  equation  relating  the  time  derivatives  of  the  stresses,  we 
add  the  equations  obtained  by  multiplying  Equation  (159)  by  u'  and 
Equation  (160)  by  x'  . 
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It  is  clear  that  Equation  (163)  can  also  be  derived  by  combining  the  time 
derivative  of  Equation  (146)  written  as 
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with  the  first  law  of  thermodynamics  expressed  by  Equation  (39) . 


Equations  (159)  through  (163)  provide  us  with  a  convenient  basis  for 
understanding  the  results  of  our  gas  gun  experiments.  The  equations 
governing  reactive  compressive  flow  produced  by  normal  impact  are 
obtained  as 
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by  setting  x£  -  du^/dh^  —  0  in  Equations  (159),  (160),  and  (163).  It 
then  follows  that  reactive  flow  produced  by  PII  is  governed  by  Equations 
(165)  and  (166)  in  regions  where  x^  -  3u0/8h^  -  0  and  by  Equations  (159) 
through  (163)  in  regions  where  x£  *  0.  We  thus  expect  the  sets  of 
longitudinal  particle  velocity  histories  recorded  in  equivalent  normal - 
impact  and  PII  experiments  to  exhibit  essentially  the  same  flow  features 
until  reaction  occurs  in  regions  influenced  by  shear  deformation. 

Further  consideration  of  these  reactive  flow  equations  leads  us  to 
conclude  that  the  flows  produced  by  normal- impact  and  PII  depend  on  the 
relative  time  scales  for  CIR,  SIR,  shear  wave  propagation,  and  the  one- 
dimensional  flow  in  our  experiments. 

Three  types  of  flow  determined  by  the  relative  magnitudes  of  these 
time  scales  are  postulated  to  provide  a  means  of  interpreting  the  results 
of  our  gas  gun  experiments.  We  label  these  flows  as  types  I,  II,  and  III 
and  characterize  them  as  follows: 

•  In  type  I  flow,  the  one -dimensional  time  of  our  gas  gun 
experiment  is  shorter  than  the  time  scales  for  CIR  and 
SIR. 

•  In  type  II  flow,  the  one -dimensional  time  of  our  gas  gun 
experiment  is  shorter  than  the  time  scale  for  CIR  but 
longer  than  that  for  SIR. 

•  In  type  III  flow,  the  one -dimensional  time  of  our  gas  gun 
experiment  is  longer  than  the  time  scale  for  CIR. 

In  the  next  two  subsections,  we  will  consider  the  occurrence  of 
these  types  of  flows  in  a  series  of  gas  gun  experiments  performed  with 
increasing  flyer  velocities  to  provide  a  framework  for  interpreting  the 
sets  of  longitudinal  particle  velocity  histories  we  have  recorded  in 
Composition  B3. 

3.3.1  Flows  Influenced  by  Compression-Induced  Reaction 

We  now  present  a  more  detailed  account  of  the  influence  of  CIR  on 
the  shock- induced  flow  produced  by  normal- impact  and  PII.  Because 
reaction  at  the  shock  front  influences  the  values  of  the  longitudinal 
particle  velocity  (u^)  and  the  compressional  stress  (P]_]_)h  there,  we 
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first  write  the  equations  for  the  time  rate  of  change  of  these  quantities 
as 
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It  follows  from  Equations  (167)  and  (168)  that  (u^)^  and  (Pq]_)y  will 

remain  unchanged  as  the  CW  propagates  if  either  (3u../3t)  -  (3u  /ah. )„  =  0 

1  H  1  I  H 

and  (ap1;L/at)H  -  (ap11/ah1)H  -  0  or  (au1/at)H  -  -  C1(au1/ah1)H  and 
(aPll/at)H  -  -  c1(ap11/ah1)H. 

Because  shock- induced  reactive  flow  is  complicated,  we  will  now 
concentrate  largely  on  a  type  of  compressive  flow  that  we  have  observed  in 
our  present  investigation  of  Composition  B3 .  In  this  type  of  flow,  which 
has  also  been  observed  in  multiple  Lagrange  gage  studies  of  other 
explosives  such  as  TNT^  and  PBX9404,^  the  shock  propagates  essentially 
unchanged  until  it  is  overtaken  by  a  reactive  stress  wave  from  the  rear. 

We  are  consequently  interested  in  shock- induced  flows  with  dpq]/3t  >  0, 
apu/ahi  <  0,  and  3u^/3t  >  0.  Although  such  flows  can  develop  initially 
with  either  3u^/3h^  <  0  or  3u^/3h^  >  0,  we  will  only  consider  the  latter 
case  here . 

It  is  convenient  to  use  the  following  equations 
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in  our  discussion  of  a  series  of  gas  gun  experiments  performed  with 
increasing  flyer  velocities.  To  provide  a  better  appreciation  of  flow 
influenced  by  CTR,  we  first  chose  the  l  ower  flyer  velocities  in  normal - 
impact  and  PI1  experiments  to  produce  flows  of  type  I .  In  this  case ,  flows 
are  nonreactive  in  the  time  frame  of  the  experiments  and  the  CWs  propagate 
as  step  shocks  with  Su^/dt  -  3u^/3h^  -  0. 

We  thus  expect  the  sets  of  longitudinal  particle  velocity  histories 
produced  at  the  lower  flyer  velocities  in  our  gas  gun  experiments  to 
exhibit  the  same  features. 

As  the  flyer  velocities  are  increased,  the  shock  becomes  strong  enough 

to  initiate  CIR  in  the  time  frame  of  the  experiment  and  we  must  consider 

flows  of  type  III.  Remembering  that  we  are  primarily  interested  in  type 

III  flow  with  the  energy  release  rate  essentially  zero  at  the  shock  front 

and  3u^/3h^  >0,  we  first  consider  the  flows  in  the  CWs  where  x£  -  du^/dh^  - 

3p2j./dt  ”  0-  Because  3u^/3h^  -  0  initially,  it  follows  from  Equation 

(159)'  that  the  stress  increases  along  a  particle  path,  3p-^^/3t  >0,  as 

the  reaction  rate  increases.  Thereafter,  because  3u^/3h^  >  0,  the  stress 

continues  to  increase  along  a  particle  until  it  attains  a  maximum, 

3pn/3t  -  0,  at  the  time  when  the  chemical  energy  release  rate  term  is 

balanced  by  the  particle  velocity  gradient  term.  Thereafter,  the  stress 

falls  along  a  particle  path  until  the  chemical  energy  release  rate 

becomes  zero  and  the  flow  becomes  isentropic.  We  next  use  Equation  (169) 

to  determine  conditions  for  u^  to  increase  along  a  particle  path  as  p^ 

increases.  Because  (3p11/3h1)„  =  0,  it  follows  from  the  momentum 

i.  1  In 

equation  and  Equation  (169)  that  u^,  will  increase  with  p^  along  a 

2 

particle  path  when  3  p^/3t3h^  <0.  We  thus  expect  u-^  to  increase  along 
a  particle  path  when  the  energy  release  rate  on  the  particle  path 
immediately  to  its  left  is  greater  than  the  energy  release  rate  on  the 
particle  path  immediately  to  its  right. 

Having  considered  the  compressive  reactive  flow  produced  in 
equivalent  normal- impact  and  PII  experiments,  we  must  now  consider  the 
additional  influence  of  CIR  on  the  flow  produced  by  PII.  Now  because 
(3p2]/<?A)e,x^,X2  *  0  when  x£  *  0,  CIR  influences  material  subjected  to 
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shear  deformation,  and  it  follows  from  Equations  (160)'  and  (163)  that 
the  energy  released  by  chemical  reaction  influences  both  3p2]_/3t  and 
3pn/3t  and  is  consequently  converted  into  the  shear  particle  velocity 
field  as  well  as  the  longitudinal  particle  velocity  field.  In  type  III 
flows  influenced  by  CIR,  we  thus  expect  more  chemical  energy  to  be 
converted  into  the  longitudinal  particle  velocity  field  in  normal- impact 
experiments  than  in  PII  experiments. 

For  the  reasons  given  in  this  discussion,  we  expect  the  sets  of 
longitudinal  particle  velocity  histories  recorded  in  equivalent  normal - 
impact  and  PII  experiments  to  exhibit  the  same  type  III  flow  features 
until  the  material  subjected  to  shear  strain  in  the  PII  experiment  is 
influenced  by  chemical  reaction.  In  the  absence  of  SIR,  differences 
between  these  sets  of  records  occur  because  chemical  energy  is  parti¬ 
tioned  between  the  longitudinal  and  shear  particle  velocity  fields  in  the 
PII  experiment  but  not  in  the  normal -impact  experiment.  We  thus  conclude 
that  qualitative  differences  between  sets  of  longitudinal  particle 
velocity  histories  recorded  in  equivalent  normal- impact  and  PII 
experiments  that  exhibit  type  III  flow  cannot  be  taken  as  evidence  for 
the  occurrence  of  SIR. 

The  influence  of  SIR  on  type  III  flows  will  be  discussed  in  the  next 
subsection. 

We  finally  note  that  type  III  flows  exhibiting  different  features 
will  be  produced  as  the  flyer  velocity  is  increased,  because  the  shock 
becomes  significantly  influenced  by  the  energy  release  rate  at  the  shock 
front . 

3.3.2  Flows  Influenced  by  Shear -Induced  Reaction 

From  our  classification  of  flows  produced  in  our  gas  gun  experiments, 

it  is  clear  that  SIR  will  first  occur  in  flows  of  type  II  produced  by 

fliers  with  velocities  too  low  to  initiate  CIR.  We  thus  consider  flows 

exhibiting  a  nonreactive  CW,  (3p../3A)  ,  ,  3A/3t  -  0,  and  chemical 

il  e ,x^ ,X£ 

reaction  initiated  by  the  shear  deformation  resulting  from  PII.  In  this 
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follows  from  Equations  (160)'  and  (162)  that  the  shear  stress  starts  to 
increase  along  a  particle  path,  3p2]_/3t  >  0.  Moreover,  because  3u^/3h^ 
is  initially  zero  but  u1  #  0  in  the  region  where  x^  +  0 ,  it  follows  from 
Equations  (159)'  and  (161)  that  the  normal  stress  p^  also  starts  to 
increase  along  a  particle  path,  dp^/dt  >  0.  These  changes  in  P2^  and 
p-^,  produced  along  a  particle  path  by  the  energy  liberated  by  SIR,  are 
related  by  Equation  (163).  It  then  follows  that  the  stress  gradients 
3pn/3hi  and  dp2\/^\,  resulting  from  the  spatial  variations  of  3p-^/3t 
and  3p2^/3t,  will  produce  corresponding  changes  in  u-^  and  U2  governed  by 
the  momentum  equations.  We  thus  expect  the  sets  of  longitudinal  particle 
velocity  histories  recorded  in  such  equivalent  normal- impact  and  PII 
experiments  to  exhibit  the  same  features  until  the  flow  produced  by  PII 
is  influenced  by  SIR.  We  believe  that  we  have  just  formulated  a  more 

3 

precise  but  qualitative  account  of  Gupta's  approach  to  the  problem  of 
SIR  in  explosives. 

We  now  return  to  type  III  flows  produced  in  normal- impact  and  PII 
experiments  and  consider  the  case  when  the  reaction  rate  is  enhanced  by 
the  shear  deformation.  It  follows  from  our  discussion  of  Equation  (136) 
that  such  enhancement  is  possible  in  our  reactive  solid  when  the  reaction 
rate  is  primarily  temperature  dependent,  because  most  of  the  shear  work 
is  converted  into  temperature.  It  also  follows  from  our  previous 
discussion  in  this  section  that  the  additional  chemical  energy  released 
in  the  PII  experiment  will  influence  both  the  longitudinal  and  shear 
particle  velocity  histories.  Thus,  we  again  expect  the  sets  of 
longitudinal  particle  velocity  histories  exhibiting  type  III  flow  to 
differ,  because  chemical  energy  is  partitioned  between  the  longitudinal 
and  shear  particle  velocity  fields  in  the  PII  experiment  but  not  in  the 
normal  impact  experiment. 

We  thus  conclude  that  evidence  for  SIR  cannot  readily  be  obtained  by 
comparing  sets  of  longitudinal  particle  velocity  histories  from  normal - 
impact  and  PII  experiments  that  exhibit  type  III  flow,  because  differences 
between  these  sets  of  histories  can  be  attributed  to  the  influence  of 
either  CIR  or  SIR  on  the  flow  produced  by  PII. 
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4.  EXPERIMENTAL  METHODS  AND  RESULTS 


This  section  describes  the  methods  used  to  prepare  the  Composition 
B3  targets  and  the  results  of  the  normal  and  parallel- inclined  impact 
experiments . 

4.1  MULTIPLE  LAGRANGE  GAGE  COMPOSITION  B3  TARGETS 

Targets  of  Composition  B3  (60%  RDX,  40%  TNT  by  weight)  were  used  in 
all  the  gas  gun  experiments.  These  targets  were  made  as  pucks  by  cast¬ 
ing  Composition  B3  around  sets  of  Lagrange  particle  velocity  gages  held 
in  a  gage  assembly.  The  casting  was  performed  at  SRI  International's 
Test  Facility,  and  preliminary  castings  were  sectioned  to  verify  that 
the  casting  procedure  produced  uniform  Composition  B3  with  no  cavities 
and  with  the  gages  firmly  embedded  in  the  explosive.  The  gage 
assemblies  used  for  casting  four  cylindrical  targets  40  mm  in  diameter 
and  20  mm  high,  each  containing  three  Lagrange  gages,  were  made  at  SRI. 
Figure  4-1  is  a  photograph  of  a  gage  assembly  showing  three  particle 
velocity  gages  before  casting.  These  gages  were  positioned  0.5,  2.5, 
and  4.5  mm  above  the  base  plate  of  the  assembly,  which  interfaces  with 
the  impact  surface  of  the  explosive  target.  The  10-mm-long  active  gage 
elements  were  made  from  0.05 -mm- diameter  copper  wires  and  carefully 
soldered  to  0 . 25-mm-diameter  copper  posts.  After  the  casting  process, 
the  copper  posts  came  out  of  the  back  surface  of  the  explosive  target, 
as  indicated  in  Figure  4-1.  The  explosive  targets  were  backed  with 

6-8  mm  of  epoxy  before  they  were  potted  in  the  Lucite  target  ring 

required  for  the  gas  gun  experiment. 

4.2  GAS  GUN  EXPERIMENTS 

We  performed  eighteen  gas  gun  experiments  with  five  batches  of 
Composition  B3.  These  shots  were  peformed  in  sets  and  were  labelled  in 

chronological  order  as  follows:  (84-2-26,  84-2-27,  84-2-28),  (84-2-33, 
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RP-6834-1 

Figure  4-1 .  Photograph  of  gage  assembly  showing  three  particle  velocity 
gages  before  the  casting  process. 
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84-2-34,  84-2-35,  84-2-36),  (85-2-8,  85-2-9,  85-2-10,  85-2-11),  (86-2- 
16,  86-2-17,  86-2-18),  and  (86-2-64,  86-2-65,  a6-2-66,  86-2-67).  The 
results  of  these  experiments  are  summarized  in  Table  4-1. 

We  will  discuss  these  sets  of  experiments  in  chronological  order 
but  will  not  present  all  the  particle  velocity  histories  because  some 
experiments  are  repeats  of  others  and  because  some  of  the  particle 
velocity  histories  are  too  noisy  to  be  useful. 

4.2.1  The  First  Set  of  Shots 

The  first  set  of  experiments  was  designed  to  produce  an  initial 
compressional  stress  in  Composition  B3  in  the  20-kbar  region.  The 
magnetic  field  was  oriented  to  measure  the  longitudinal  particle 
velocity  histories  produced  by  normal  impact  in  Shot  84-2-26  and  by  PII 
(15°  angle)  in  Shot  84-2-28  and  the  shear  particle  velocity  histories 
produced  by  PII  (15°  angle)  in  Shot  84-2-27.  Figure  4-2(a)  shows  a 
composite  plot  of  the  longitudinal  particle  velocity  histories  produced 
in  Shot  84-2-26  by  the  normal  impact  of  a  1/4- inch- thick  fused  silica 
flyer  travelling  at  0.558  mm/ns ,  and  Figure  4- 2(b)  shows  a  composite 
plot  of  those  produced  in  Shot  84-2-28  by  the  PII  of  a  3/8- inch- thick 
fused  silica  flyer  travelling  at  0.575  mm/ ns.  When  interpreting  these 
particle  velocity  histories,  it  is  important  to  remember  that  the 
experimental  frame  time  at  a  gage  location  is  -3.5  /is .  Data  beyond 
3.5  p is  are  questionable  because  of  lead  effects. 

We  now  recall  that  our  investigation  of  SIR  is  based  on  a 
comparison  of  the  sets  of  longitudinal  particle  velocity  histories 
produced  by  similar  compressional  loading  conditions  in  normal- impact 
and  PII  experiments.  Sets  of  histories  were  compared  by  superposition¬ 
ing  one  set  of  histories  on  the  other.  Superposition  of  the  histories 
in  Figures  4-2(a)  and  4-2(b)  showed  (1)  that  these  sets  of  histories 
cannot  be  compared  at  all  three  gage  positions  because  the  CW  broke  the 
third  gage  in  the  normal- impact  experiment,  and  (2)  that  the  historie: 
recorded  by  the  first  and  second  gages  in  both  types  of  impact  experi¬ 
ment  are  essentially  the  same  when  the  flyer  plate  thicknesses  are  taken 
into  account.  Comparison  of  the  histories  in  Figure  4-2 (b)  reveals 
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Table  4-1 


SUMMARY  OF  COMPRESSION- SHEAR  EXPERIMENTS 
PERFORMED  IN  COMPOSITION  B3a 


Shot 

No. 

Flyer  Velocity 
(mm/MS ) 

Angle 

Degrees 

Initial  Particle 
Velocity  (mm//is) 

Comments 

84-2-26b 

0.558 

0° 

0.416 

Third  gage  broke 

84-2-27c 

0.579 

15° 

None 

No  good  records 

84-2  - 28b 

0.575 

15° 

0.414 

Reactive  flow 

84-2-33b 

0.295 

0° 

0.215 

Reactive  flow 

84-2-34b 

0.292 

15° 

0.220 

Reactive  flow 

84-2-35° 

0.302 

15° 

None 

No  good  records 

84- 2  - 36b 

0.300 

IS¬ 

None 

First  gage  record 
too  noisy 

85-2-8b 

0.301 

C'1 

0.216 

Reactive  flow 

85-2-9b 

0.305 

15° 

0.200 

Reactive  flow 

85-2-10b 

0.308 

20° 

0.200 

Reactive  flow 

85-2-11° 

0.312 

15° 

None 

No  good  records 

86-2-16b 

0.294 

0° 

None 

First  gage  record 
too  noisy 

86-2-17b 

0.304 

15* 

0.140 

Essentially 

nonreactive 

86-2-18b 

0.294 

15° 

0.136 

Second  gage  faulty 

86  -  2  - 64b 

0.293 

0° 

0.230 

Essentially 

nonreactive 

86-2-65b 

0.302 

15° 

0.230 

Essentially 

nonreactive 

86-2-66b 

0.452 

15° 

0.310 

Reactive  flow 

86-2-67b 

0.429 

0° 

0.320 

Reactive  flow 

aA  1/4- inch- Chick  fused  silica  flyer  was  used  in  Shot  84-2-26;  3/8  -  inch- thick 
fused  silica  flyers  were  used  in  all  the  others. 

^Magnetic  field  aligned  to  measure  longitudinal  particle  velocity. 
cMagnetic  field  aligned  to  measure  shear  particle  velocity. 


50 


PARTICLE  VELOCITY  (mm/ps)  PARTICLE  VELOCITY  (mm/ps) 


0.0  2.0  4.0  6  0 

TIME  (ps) 

(a)  Shot  84-2-26,  normal  impact 


0.0  2.0  4.0  6  0 

TIME  (ps) 

(b)  Shot  84-2-28,  parallel-inclined  impact  (15°  angle) 

J  A-6S3J-1 

Figure  4-2.  Composite  plots  of  Lagrange  longitudinal  particle  velocity 
histories  in  Composition  B3  produced  by  normal  impact  in 
Shot  84-2-26  and  by  parallel-inclined  impact  (15°  angle)  in 
Shot  84-2-28. 
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developing  flow  features  characteristic  of  a  type  III  flow,  produced  by 

a  shock  strong  enough  to  initiate  a  significant  energy  release  rate  at 

the  shock  front.  The  particle  velocity  (u-^)^  at  the  shock  front 

increases  d(u^)^/dt  >  0,  and  until  the  flow  behind  the  shock  is 

influenced  by  the  rarefaction  from  the  rear  of  the  flyer  plate,  the 

compressional  stress  gradient,  dp^/^h^,  is  zero  at  the  first  and  second 

gage  positions,  where  3u^/3t  -  0,  but  is  negative  at  the  third  gage 

position  where  d\i-/8t  >  0.  Because  d(u- )  /dt  >0,  it  follows  from 

Equation  (166)  that  the  energy  release  rate  at  the  shock  front  is 

sufficient  to  make  d(p^-^)^/dt  >  0  and  (dp^/dt)^  >0.  In  addition,  the 

energy  released  by  CIR  behind  the  shock  front  causes  p^  to  increase 

along  a  particle  path  in  such  a  way  that  the  second  derivative 
2 

(3  p^/3t3h.)  is  zero  in  the  neighborhood  of  the  first  two  gages  but 
negative  in  the  neighborhood  of  the  third  gage. 

We  conclude  from  this  discussion  that  our  first  set  of  shots  is 
inconclusive  with  regard  to  SIR  because  the  initial  compressional 
loading  in  the  20-kbar  region  produces  a  type  III  flow  that  is 
significantly  influenced  by  CIR. 

The  records  from  Shot  84-2-27,  the  experiment  designed  to  produce 
shear  particle  velocity  histories,  are  not  presented  because  no 
significant  shear  particle  velocity  was  recorded  beyond  the  first  gage 
position. 

4.2.2  The  Second  Set  of  Shots 

The  second  set  of  experiments  was  designed  to  produce  an  initial 
compressional  stress  in  Composition  B3  in  the  10-kbar  region.  The 
magnetic  field  was  oriented  to  measure  the  longitudinal  particle 
velocity  histories  produced  by  normal  impact  in  Shot  84-2-33  and  by 
PII  (15°  angle)  in  Shots  84-2-34  and  84-2-36  and  the  shear  particle 
velocity  histories  produced  by  PII  (15°  angle)  in  Shot  84-2-35.  Shot 
84-2-36  was  performed  as  a  repeat  of  Shot  84-2-34  to  test  the  repro¬ 
ducibility  of  our  experimental  technique,  but  the  excessive  noise 
recorded  at  some  of  its  gage  positions  invalidated  this  test.  This  is 
also  the  reason  why  we  do  not  present  the  records  from  Shot  84-2-36. 
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A  3/8 -inch- thick  fused  silica  flyer  was  used  in  all  of  these  shots. 
Figure  4-3(a)  shows  a  composite  plot  of  the  longitudinal  particle 
velocity  histories  produced  in  Shot  84-2-33  by  the  normal- impact  of  a 
flyer  travelling  at  0.295  mm/^s ,  and  Figure  4- 3(b)  shows  a  composite 
plot  of  those  produced  in  Shot  84-2-34  by  the  PII  (15°  angle)  of  a  flyer 
travelling  at  0.292  mm/jus.  The  voltage  traces  recorded  in  Shot  84-2-34 
by  the  analog  oscilloscopes  are  also  presented  (Figure  4-4)  to  exemplify 
the  quality  of  the  particle  velocity  data  usually  recorded  in  our 
experiments . 

Superposition  of  the  histories  in  Figures  4-3(a)  and  4-3(b)  shows 
(1)  that  the  sets  of  histories  can  be  compared  at  all  three  gage 
positions,  and  (2)  that  the  flows  recorded  in  Shots  84-2-33  and  84-2-34 
initially  developed  in  the  same  way  but  then  became  significantly 
different.  Properties  of  the  particle  velocity  histories  recorded  in 
both  shots  are  discussed  in  some  detail  below  to  specify  the 
similarities  and  differences  between  the  two  flows. 

Qualitatively,  the  particle  velocity  histories  recorded  at  the 
first  gage  plane  are  similar  and  exhibit,  behind  the  wave  front,  a 
region  of  constant  particle  velocity  followed  by  a  region  where  the 
particle  velocity  rises  to  a  peak  and  then  falls.  Quantitatively,  these 
histories  are  essentially  identical  for  1.76  ns  and  then  become 
different.  After  1.76  ns ,  the  particle  velocity  recorded  in  Shot 
84-2-34  rises  to  peak  value  of  0.275  mm//is  at  3.0  ns,  but  the  particle 
velocity  recorded  in  Shot  84-2-33  remains  essentially  constant  for 
another  0.5  ns  before  rising  to  a  peak  value  of  0.239  mm//is  at  2.78  ns  ■ 

Qualitatively,  the  particle  velocity  histories  recorded  at  the 
second  gage  plane  are  different  in  the  sense  that  the  history  from  Shot 
84-2-34  exhibits  a  well-defined  particle  velocity  peak  but  the  history 
from  Shot  84-2-33  does  not.  Quantitatively,  however,  these  histories 
are  essentially  identical  for  1.06  ns ■  After  1.06  ns,  the  particle 
velocity  recorded  in  Shot  84-2-34  rises  to  a  peak  value  of  0.219  mm/ ns 
at  2.12  ns ■  but  the  particle  velocity  recorded  in  Shot  84-2-33  remains 
essentially  constant  for  most  of  the  history. 
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Figure  4-4.  Electromagnetic  particle  velocity  gage 
records  from  the  analog  oscilloscopes 
in  Shot  84-2-34. 
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The  particle  velocity  histories  recorded  at  the  third  gage  plane 
are  both  qualitatively  and  quantitatively  different.  The  particle 
velocity  in  Shot  84-2-34  rises  to  a  value  of  0.2  mm//us  at  the  wave  front 
and  remains  constant  for  2  /us.  However,  the  particle  velocity  recorded 
in  Shot  84-2-33  rises  to  a  value  of  0.207  mm//us  at  the  wave  front, 
continues  to  rise  and  attains  a  value  of  0.213  mm//ts  at  2.25  /us ,  then 
falls  to  a  value  of  0.207  mm//us  at  2.68  /us,  and  finally  rises  to  a  value 
of  0.25  mm/jus  at  4.33  /us. 

Further  examination  of  these  sets  of  particle  velocity  histories 
leads  us  to  consider  two  alternative  interpretations  of  these  flows.  In 
the  first  interpretation,  we  take  the  records  at  face  value  and  see  that 
both  flows  develop  as  the  same  decaying  step  shock  in  the  (p-Q  -  h-^) 
plane  until  the  particle  velocity  starts  to  increase  in  the  PII 
experiment.  Therefore,  the  stress  gradient  becomes  negative  at  the 
first  and  second  gage  planes  in  Shot  84-2-34  and  at  the  first  gage  plane 
in  Shot  84-2-33.  We  now  use  Equation  (168)  to  show  that  (3p^^/3t)H  <  0 
at  the  front  of  the  decaying  shock  because  d(p^-^)^/dt  <  0  and 
(3pn/3hi)H  -  0,  and  to  infer  that  Sp-^/dt  <  0  and  S^p-^/StSh^  -  0  in 
the  decaying  step  shock  where  3u^/3t  -  0. 

We  can  then,  by  extending  our  previous  discussions  of  reactive 
flow,  associate  the  decaying  step  shock  with  a  negative,  global,  energy 
release  rate  and  the  ensuing  negative  stress  gradients  with  a  positive, 
global,  energy  release  rate.  This  being  the  case,  it  is  important  to 
note  (1)  that  the  negative  stress  gradient  appears  first  in  the  flow 
produced  by  PII  in  Shot  84-2-34,  and  (2)  that  the  particle  velocity  in 
Shot  84-2-34  increases  for  a  longer  period  of  time  and  attains  a 
significantly  higher  peak  value  than  that  in  Shot  84-2-33. 

In  the  second  interpretation  of  these  two  flows,  the  differences 
between  the  shock  particle  velocities  (u^)^  recorded  at  the  first  and 
second  gage  planes  in  both  shots  are  neglected  for  practical  purposes. 
Then  the  global  energy  release  rate  in  both  flows  is  considered  to  be 
initially  zero,  but  our  previous  argument  associating  the  negative 
stress  gradients  with  a  positive  energy  release  rate  is  still  valid. 
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From  these  two  interpretations  of  the  flows  recorded  in  Shots 

84- 2-33  and  84-2-34,  we  conclude  that  the  influence  of  the  energy 
release  rate  occurs  first  in  the  flow  produced  by  PII,  and  consequently, 
we  designate  this  flow  as  a  kind  of  type  II  flow.  We  then  accordingly 
attribute  the  significant  differences  between  the  sets  of  longitudinal 
particle  velocity  histories  recorded  in  Shots  84-2-33  and  84-2-34  to  the 
occurrence  of  SIR. 

The  records  from  Shot  84-2-35  are  not  presented  because  they  are 
too  noisy. 

4.2.3  The  Third  Set  of  Shots 

The  third  set  of  experiments  was  also  designed  to  produce  an 
initial  compressional  stress  in  Composition  B3  in  the  10-kbar  region. 

The  magnetic  field  was  oriented  to  measure  the  longitudinal  particle 
velocity  histories  produced  by  normal  impact  in  Shot  85-2-8,  by  PII  (15° 
angle)  in  Shot  85-2-9,  and  by  PII  (20°  angle)  in  Shot  85-2-10;  and  the 
shear  particle  velocity  histories  produced  by  PII  (15°  angle)  in  Shot 

85- 2-11.  A  3/8 -inch- thick  fused  silica  flyer  was  used  in  all  of  these 
experiments.  Composite  plots  of  the  longitudinal  particle  velocities 
recorded  in  Shot  85-2-8  (0°  angle)  with  the  flyer  travelling  at  0.301 
mm//xs ;  in  Shot  85-2-9  (15°  angle),  with  the  flyer  travelling  at  0.305 
ram/)j,s ;  and  in  Shot  85-2-10  (20°  angle),  with  the  flyer  travelling  at 
0.308  mm//is ;  are  shown  in  Figures  4-5  through  4-7.  The  results  of  the 
shear  particle  velocity  measurements  obtained  in  Shot  85-2-11  (15° 
angle),  with  the  flyer  travelling  at  0.312  mm//is ,  are  shown  in  Figure 
4-8.  These  records  from  Shot  85-2-11  are  presented  here  to  typify  the 
results  of  our  shear  particle  velocity  measurements  in  this  work  and 
demonstrate  our  inability  to  record  well-defined  shear  particle  velocity 
histories  in  Composition  B3. 

Superposition  of  the  sets  of  longitudinal  particle  velocity 
histories  in  Figures  4-5  through  4-7  leads  to  three  important  observations . 
The  first  is  that  the  flows  produced  by  normal-impact  and  PII  in  Shots 
85-2-8,  85-2-9,  and  85-2-10  are  for  practical  purposes  the  same  until 
they  are  significantly  influenced  by  the  energy  release  rate.  The  second 
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Figure  4-6.  Composite  plot  of  Lagrange  longitudinal  particle  velocity  histories  in 
Composition  B3  produced  by  parallel-inclined  impact  (15°  angle) 
in  Shot  85-2-9. 
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Figure  4-7.  Composite  plot  of  Lagrange  longitudinal  particle  velocity  histories 
in  Composition  B3  produced  by  parallel-inclined  impact  (20°  angle) 
in  Shot  85-2-10. 
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Figure  4-8.  Records  from  shot  85-2-11,  a  parallel-inclined  (15°  angle)  impact 
experiment  designed  to  measure  Lagrange  shear  particle  velocity 
histories  in  Composition  B3. 
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is  chat  the  energy  release  rate  has  a  more  pronounced  on  the  flow 
produced  in  the  normal  impact  experiment  than  on  the  flows  produced  in 
the  PII  experiments.  The  third  is  that  all  the  sets  of  longitudinal 
particle  velocity  histories  exhibit  characteristics  of  type  III  flow 
influenced  by  CIR. 

We  conclude  as  a  result  of  these  observations  that  differences  among 
the  sets  of  longitudinal  particle  velocities  recorded  in  Shots  85-2-8, 
85-2-9,  and  85-2-10  cannot  be  taken  as  evidence  for  the  occurrence  of 
SIR. 

Because  different  conclusions  were  reached  from  the  results  of  the 
second  and  third  sets  of  experiments,  and  both  sets  of  experiments  were 
designed  to  produce  initial  compressional  stresses  in  the  10-kbar  region, 
we  now  asked  why  the  sets  of  longitudinal  particle  velocity  histories 
provide  conflicting  evidence  for  the  occurrence  of  SIR.  The  only  obvious 
answer  is  that  the  sets  of  longitudinal  particle  velocity  histories 
recorded  in  the  second  and  third  sets  of  gas  gun  experiments  exhibit 
different  flow  features  because  these  sets  of  experiments  were  performed 
with  different  batches  of  Composition  B3.  This  answer  led  us  to  suspect 
that  the  batches  of  Composition  B3 ,  seemingly  cast  under  the  same 
conditions,  had  different  physical  characteristics. 

4.2.4  The  Fourth  Set  of  Shots 

Because  different  results  were  obtained  in  the  second  and  third  sets 
of  experiments,  the  fourth  set  was  again  designed  to  produce  an  initial 
compressional  stress  in  Composition  B3  in  the  10-kbar  region.  The 
magnetic  field  was  oriented  to  measure  the  longitudinal  particle 
velocities  produced  by  normal  impact  in  Shot  86-2-16  and  by  PII  (15° 
angle)  in  Shots  86-2-17  and  86-2-18.  These  experiments  were  again 
performed  with  a  3/8  -  inch- thick  fused  silica  flyer.  The  velocity  of  the 
flyer  was  0.294  mm//is  in  Shot  86-2-16,  0.304  mm/jus  in  Shot  86-2-17,  and 
0.294  mm/^s  in  Shot  86-2-18.  Figure  4-9  shows  a  composite  plot  of  the 
longitudinal  particle  velocity  histories  recorded  in  Shot  86-2-17.  The 
records  from  Shot  86-2-16  are  not  presented  because  they  are  excessively 
noisy,  and  those  from  Shot  86-2-18  are  not  presented  because  the  second 
particle  velocity  gage  was  found  to  be  inoperative  in  this  experiment. 
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Figure  4-9.  Composite  plot  of  Lagrange  longitudinal  particle  velocities 

in  Composition  B3  produced  by  parallel-inclined  impact  (15:  angle) 
in  Shot  86-2-17. 


Examination  of  the  composite  plot  in  Figure  4-9  showed  that  a  flow 
of  type  I  exhibiting  little  evidence  of  CIR  and  no  evidence  of  SIR  is 
produced  in  Shot  86-2-17.  Comparison  of  the  particle  velocity  histories 
in  Figure  4-9  with  those  in  Figures  4-3(b)  and  4-6,  which  were  produced 
by  similar  flyer  velocities,  then  showed  that  anomalously  low  particle 
velocities  were  recorded  in  Shot  86-2-17.  This  anomalous  result  and  the 
observation  from  Figure  4-9  that  the  formation  of  a  precursor  wave  is 
more  evident  in  Shot  86-2-17  than  in  any  of  the  other  experiments  were 
taken  as  further  evidence  that  Composition  B3  targets  prepared  in 
different  castings  had  different  physical  properties. 

4.2.5  The  Fifth  Set  of  Shots 

Because  anomalous  particle  velocities  were  recorded  in  the  fourth 
set  of  experiments  and  flows  of  type  III  were  produced  in  the  first  set 
of  experiments,  the  fifth  set  of  experiments  was  designed  to  produce 
initial  compressional  stresses  in  Composition  B3  in  the  10-  and  15-kbar 
regions.  The  magnetic  field  was  oriented  to  measure  the  longitudinal 
particle  velocities  produced  by  normal  impact  in  Shots  86-2-64  and 
86-2-67  and  those  produced  by  PII  (15°  angle)  in  Shots  86-2-65  and 
86-2-66.  A  3/8 -  inch- thick  fused  silica  flyer  was  again  used  in  all  of 
these  experiments.  Composite  plots  of  the  longitudinal  particle 
velocities  recorded  in  Shot  86-2-64  (0°  angle),  with  the  flyer  travelling 
at  0.293  mm//is ;  in  Shot  86-2-65  (15°  angle),  with  the  flyer  travelling  at 
0.302  mm/ns ;  in  Shot  86-2-66  (15°  angle)  with  the  flyer  travelling  at 
0.452  mm/zis ;  and  in  Shot  86-2-67  (0°  angle),  with  the  flyer  travelling  at 
0.429  mm/ps,  are  shown  in  Figures  4-10  through  4-13. 

Superposition  of  the  composite  plots  in  Figures  4-10  and  4-11  shows 
that  the  particle  velocity  histories  recorded  in  Shots  86-2-64  and 
86-2-65  are  essentially  identical  and  show  little  evidence  of  exothermic 
reaction.  Further  examination  of  these  sets  of  particle  velocity 
histories  again  leads  to  two  alternative  interpretations  of  these  flows; 
these  interpretations  have  already  been  discussed  (see  Section  4.2.2)  in 
our  consideration  of  the  particle  velocity  histories  shown  in  Figures 
4-3(a)  and  4-3(b).  Here  for  pract.  xl  purposes  we  neglect  the  decaying 
step  shock  associated  with  a  negative  energy  release  rate  and  claim  that 
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Figure  4-10.  Composite  plot  of  Lagrange  longitudinal  particle  velocity  histories 
in  Composition  B3  produced  by  normal  impact  in  Shot  86-2-64. 
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Figure  4-1 1 .  Composite  plot  of  Lagrange  longitudinal  particle  velocity  histories 
in  Composition  B3  produced  by  parallel-inclined  impact  (15:  angle) 
Shot  86-2-65. 
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Figure  4-12.  Composite  plot  of  Lagrange  longitudinal  particle  velocity  histories 
in  Composition  B3  produced  by  parallel-inclined  impact  (15;  angle) 
in  Shot  86-2-66. 
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Figure  4-13.  Composite  plot  of  Lagrange  longitudianl  particle  velocity  histories 
in  Composition  83  produced  by  normal  impact  in  Shot  86-2-67. 
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the  flows  produced  by  normal- impact  and  PII  in  Shots  86-2-64  and  86-2-65 
are  nonreactive  flows  of  type  1. 

Superposition  of  the  composite  plots  in  Figures  4-12  and  4-13  shows 
that  the  flows  produced  by  normal- impact  and  PII  in  Shots  86-2-67  and  86- 
2-66  are  initially  identical  (for  about  2.3  /is)  ,  but  then  become  markedly 
different.  After  2.3  ns,  the  stress  gradient  at  the  first  gage  position 
becomes  negative  in  Shot  86-2-66  but  positive  in  Shot  86-2-67,  and  this 
difference  between  the  flows  produced  by  PII  and  normal  impact  cannot  be 
neglected.  Thereafter,  the  particle  velocity  histories  recorded  at  the 
second  and  third  gage  positions  in  both  experiments  exhibit  significant 
reactive  flow  features  typical  of  a  type  III  flow,  but  the  reaction  has  a 
more  pronounced  influence  in  the  flow  produced  by  normal  impact  than  in 
that  produced  by  PII. 

In  the  light  of  our  theoretical  discussions  about  reactive  flows 
produced  in  gas  gun  experiments ,  we  conclude  that  these  differences 
between  the  sets  of  longitudinal  particle  velocity  histories  recorded  in 
Shots  86-2-66  and  86-2-67  can  be  taken  as  evidence  for  SIR  produced  by 
PII. 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 


Experimental  and  theoretical  studies  of  compression  and  shear  wave 
propagation  in  Composition  B3  were  performed  to  establish  the  occurrence 
of  shear- induced  reaction  (SIR)  and  provide  a  better  understanding  of  its 
influence  on  the  one -dimensional  flow  produced  in  an  explosive  by  the 
type  of  parallel-inclined  impact  (PII)  experiment  pioneered  by  Gupta. ^ 

Our  theoretical  study  was  based  on  a  model  analysis  of  the  PII 
experiment  with  the  shock- loaded  target  treated  as  an  inert  elastic  solid 
and  the  resulting  flow  treated  ideally  as  a  pure  compression  wave 
followed  by  a  pure  shear  wave.  A  complete  equation  of  state  (EOS)  for  a 
reactive  elastic  solid  was  formulated  and  used  to  derive  the  thermo¬ 
dynamic  equations  governing  the  one -dimensional  adiabatic  flows  produced 
by  PII.  The  conditions  an  elastic  solid  subjected  to  PII  must  satisfy  to 
support  our  ideal  inert  flow  were  then  formulated  and  used  to  construct 
the  most  convenient  model  EOS  for  investigating  the  reactive  adiabatic 
flows  produced  in  solids  by  normal- impacts  and  PIIs.  Flows  influenced  by 
compression- induced  reaction  (CIR)  and  by  SIR  were  then  considered  to 
show  how  the  energy  released  by  chemical  reaction  influences  the 
longitudinal  particle  velocity  field.  The  results  of  this  consideration 
of  reactive  flow  in  a  solid  were  then  used  to  formulate  a  framework  for 
interpreting  the  results  of  our  gas  gun  experiments. 

Our  experimental  stuiy  was  based  on  Gupta's  premise  that  significant 
differences  between  the  sets  of  longitudinal  particle  velocity  histories 
recorded  in  an  explosive  subjected  to  normal- impact  and  PII  producing 
equivalent  nonreac.tive  compression  waves  can  be  taken  as  evidence  for  the 
occurrence  of  SIR  in  the  PII  experiment.  Eighteen  normal- impact  and  PII 
experiments  were  performed  in  a  gas  gun  with  Composition  B3  targets,  and 
multiple  Lagrange  gages  embedded  in  these  explosive  targets  were  used  to 
make  longitudinal  and  shear  particle  velocity  measurements.  The  eighteen 
experiments  were  divided  into  five  sets  that  were  performed  with  targets 
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made  from  different  batches  of  Composition  B3 .  The  results  of  the 
longitudinal  particle  velocity  measurements  gave  evidence  for  the 
occurrence  of  SIR  in  the  first  batch  of  Composition  B3,  which  was  shocked 
initially  into  the  10-kbar  region,  and  in  the  fifth  batch  of  Composition 
B3 ,  which  was  shocked  initially  into  the  15-kbar  region.  The  results  of 
the  shear  particle  measurements  led  to  the  conclusion  that  the  shear 
deformation  produced  in  Composition  B3  by  PII  remains  close  to  the  impact 
surface . 

We  took  a  simplistic  approach  to  the  propagation  of  compression  and 
shear  waves  in  explosives  by  treating  the  explosive  as  an  elastic  solid 
governed  by  an  equation-of-state  that  decouples  the  longitudinal  and 
shear  deformations  when  the  flow  is  nonreactive .  This  approach  is 
relevant  in  the  light  of  our  experimental  approach  to  SIR  because  it 
allows  us  to  concentrate  on  reactive  flow  and  the  influence  of  chemical 
reaction  on  the  longitudinal  particle  velocity  field.  It  is  clear  that 
such  an  understanding  of  reactive  flow  is  required  to  interpret  the  sets 
of  longitudinal  particle  velocity  histories  recorded  in  equivalent 
normal- impact  and  PII  experiments  and  attribute  significant  differences 
in  their  reactive  flow  features  to  SIR. 

The  results  of  our  theoretical  study  substantiate  and  provide  a  more 
detailed,  qualitative  account  of  Gupta's  experimental  approach  to  SIR, 
but  they  also  show  that  the  complex  nature  of  reactive  flow  makes  this 
approach  more  complicated  than  was  originally  thought.  A  more  quanti¬ 
tative  procedure  is  needed  for  interpreting  the  sets  of  longitudinal 
particle  velocity  histories  recorded  in  the  gas  gun  experiments. 
Additional  work  is  also  required  to  understand  why  no  well-defined  shear 
particle  velocity  histories  were  recorded  in  our  PII  experiments  with 
Composition  B3 . 

Although  the  results  of  our  gas  gun  experiments  with  Composition  B3 
provide  evidence  for  the  occurrence  of  SIR  produced  by  PII,  they  are 
inconsistent  and  therefore  not  totally  convincing.  Because  experimental 
conditions  were  essentially  the  same  in  the  second,  third,  and  fourth 
sets  of  experiments,  inconsistencies  among  the  sets  of  longitudinal 
particle  velocity  histories  recorded  in  these  experiments  must  be 


72 


attributed  to  differences  among  the  Composition  B3  targets.  It  thus 
follows  that  our  casting  procedure  produced  batches  of  Composition  B3 
targets  with  different  physical  properties  and  was  not  as  well  controlled 

,  as  we  had  originally  thought.  More  attention  must  be  given  to  the 

casting  procedures  used  to  make  targets  for  the  gas  gun  experiments. 

*  We  conclude  that  the  results  of  our  theoretical  and  experimental 
studies  provide  evidence  for  the  occurrence  of  SIR  in  explosives  subject 
to  PII,  but  not  enough  to  establish  that  SIR  plays  a  significant  role  in 
the  shock  initiation  of  explosives.  We  recommend  that  any  future 
theoretical  study  of  SIR  be  based  on  a  more  realistic  description  of  an 
explosive.  This  description  can  then  be  used  to  develop  a  more  quanti¬ 
tative  procedure  for  interpreting  the  sets  of  longitudinal  particle 
velocity  histories  recorded  in  equivalent  normal- impact  and  PII 
experiments  and  also  to  determine  why  no  well-defined  shear  particle 
velocity  histories  were  recorded  in  our  PII  experiments.  We  finally 
recommend  that  any  future  gas  gun  investigations  of  SIR  be  performed  with 
explosive  targets  with  similar  physical  properties  produced  by  a  well- 

* 

proved  casting  technique. 
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